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SUMMARY 


The control of large space structures encompasses a multitude of physical 
phenomena. The structure Itself Is a complex vibrating system that Is excited 
by internal and external forces. The external forces and torques come from 
aerodynamic, sc ; nr wind, and thermal excitations to geomagnetic and gravity 
gradient forces. Internal forces and torques are created by vibrating machinery 
(CMC 8, gyros, etc.), by articulating structural elements and motions of astronauts. 
It is not surprising, then, that as the size and performance demands on structures 
increase the control problem looms ever larger as one of the overriding problems. 

This final report describes the work that was performed by the Grumman 
Aerospace Corporation Research Department under contract to the Marshall Space 
Flight Center (Contract NAS 8-32587) which was administered by Dr. Michael 
Borelli. The thrust of this effort was to determine what, if any, limits ons are 
imposed on the size of spacecraft which may be controlled using current control 
system design technology. The particular problems investigated were: 

1. The fundamental limitations imposed by structural/control 
interactions, by external torques, and by the mission 
performance requirements for Low Earth Orbit (LEO) missions. 

2. The development of control approaches for the various control 
tasks that are required by large space structures, i.e,, 

as required during fabrication, assembly, pointing, shape 
and attitude control, etc. 

3. The development of techniques for on orbit dynamic testing that 
will permit evaluation, during operation, of the parameters re- 
quired for control design. 

4. Investigate actuator requirements so that the control may be 
achieved with minimal use of expendable fuels. 

These tasks were investigated by using a typical structure in the 35 to 70 
meter size category. A control system design that used actuators that are 
currently available (CMG's) was designed for this structure. The amount of 
control power required to maintain the vehicle in a stabilized gravity gradient 


ill 


pointing orientation that also damped various structural motion t was determined. 

The moment of Inertia and mass properties of this structure were varied to 
verify that stability and performance were maintained. This was accomplished by 
systematically varying the linear dimensions of the structure by a scale factor 
while maintaining the control gains fixed. The conclusion from this study was 
that the structure's size was required to change by at least a factor of two before 
any stability problem arose. The stability margin that was lost was due to 
the scaling of the gravity gradient torques (the rigid body control) and as 
such could easily be corrected by changing the control gains assoclattid with 
the rigid body control. A secondary conclusion from this study was that the 
control design that accommodates the structural motions (to damp them) Is a 
little more sensitive than the design that works on attitude control of the 
rigid body only. The main details of this effort are described In Section I. 

The control of large structures can be considered, as the classical control 
system designer does, as the problem of controlling the attitude of the vehicle 
despite all of the disturbances that are potentially misorlentatlng the vehicles 
or causing stability problems. This approach artlflcally divides the physical 
phenomena Into control forces and disturbance forces. VJhen this Is done, the 
possibility of balancing one disturbance by a second Is lost. Thus, gravity- 
gradient forces, which are oscillatory, might be damped by using residual 
aerodynamic forces to extract energy from the gravity gradient mode. To be 
able to achieve this goal, early In the program It was perceived that a fairly 
comprehensive model of the spacecraft was needed. This model (described In 
Section II) was used to develop a control system that uses a linear optimal 
control approach to: 

• Achieve a stable rigid body orientation using gravity gradient 
forces balanced by residual aerodynamic forces 

f Stabilize the interaction of rigid and flexlbJ';; motion 

• Increase the damping of the more Important lllexlble modes 
using the same actuators that are used for rigid body 
control (a set of three orthogonal CMG's only) 

• Be slewed, when required, from one orientation to another 
without adversely affecting the vibration and shape of the 
structure 
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As the control study developed, several problems emerged 
that were either poorly understood or that were critically lacking in theoretical 
basis. These ware: 

• A method is required for order reduction of the structural 
dynamics that accounts for the undamped vibration modes that are 
left out and also considers the effect of control bandwidth that 
exceeds the frequencies: of the modes neglected 

• A method is required that reduces the control excitation of the 
high frequency modes. 

• A method is required that permits measurements that are not 
adversely corrupted by unwanted structural oscillations 

• It became clear that an estimator is undesirable in the 
control loop because the bandwidth of the filter might easily 
dominate the problem. The resulting loss of gain and phase margin 
compared with an optimal control design is clearly undesirable. We 
instead have proposed that control designs be undertaken which 
utilize ar many rate and position measurements as there are 

modes retained 

• A method for on-orblt dynamic testing Is required 

In Section III, these problems are discussed and some solutions are proposed. 
We have shown; 

• A new method for order reduction which both Incorporates the 
closed loop dynamic characteristics and the unique problem 
associated with finite element modeling which Ignores structural 
damping 

m A method for control spillover reduction at higher frequencies 
using a low order observer 

• A method for on-orblt dynamic testing which gives structural mode 
data and also reduces the measurement spillover problem 

• A new method for synoptic control design that naturally suggests 
alternate actuators 
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I, INTRODUCTION? 


As the size of space structures Increase and the performance requirements 
become more severe the control problem looms larger and larger as the single 
most difficult problem that must be overcome. The structural size alone 
creates a possibility for interactions between the control system and the 
structural dynamics. The control performance requirements on both pointing and 
Internal vibration (Jitter) create a need for new and novel techniques for rigid 
body control that simultaneously achieves the desired pointing accuracy vd.thout 
creating structural vibration that might jitter crucial spacecraft systems 
(optics* antennas, etc.). The rigid body control must also be achlevid in 
such a way that internal strain on the spacecraft structure does not change at 
certain desired points (points where critical optical or antenna elements are 
mounted must always remain at the same relative positions so that the phasing 
of the light or radio freq&4Ucy energy does not become distorted). The achievement 
of this goal requires control systems with zero stesf y state error after a slew 
command. The structures that are being proposed ch/nge their size as a 
function of time when the structure is be.ing fabricated In orbit. The ability 
f-o alter the control system as the dynamics change, the ability to recognize 
the way in which the structural and rigid body dynamics change, and the ability 
to sense the disturbances and the way they change are all desirable characteristics 
that the control system should possess. 

To achieve the flexibility that all of the above requirements Impose on ^ he 
control system requires a fairly complex system, one which also has Inherent 
reliability and insensitivity to variations in the structural dynamics. Also 
one must attempt to exploit the existing physical phenomena that cause "disturb- 
ances’* on the spacecraft to control the spacecraft. Thus gravitational, geo- 
magnetic, residual aerodynamic and solar wind forces that are normally considered 
as disturbances that are countered by expending reaction jet fuel, should be 
designed into the control system to provide the possibility of "playing one 
force off against another". 

To achieve all of these goals a synoptic approach has been developed that 
attempts to design the control system for. the spacecraft so that, at the outset, all 
of the dominant effect are modeled - the spacecraft dynamics are modeled, using 
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a finite element approach* with both the rigid body and flexible motion coupled 
In a coordinate eyate^ that doee not necesaarlly explicitly Include the rigid 
body - the phyaical forces that interact with the structure are Included explicitly 
08 distributed forces and torquea acting on each of the masses (or Inertias) 

In the finite element model — the Internal vibration sources are modeled and 
Included to provide a measure of the Jitter Induced by these motions When they 
are deemed critical — the slew commands are included in the dynamic description 
so that the best possible method for commanding the system to change its orien- 
tation can be determined — the flexibility Inherent with digital control 
requires that all of the modeling be done in such a way that any problems 
introduced by sampling can be minimized. 

The study that is described in this report has developed the technology 
to address some of these questions and has pointed to problems where technology 
must be further developed. Figure 1 shows the structure that was used to 
develop the control technology and evalu 3 ?,te the effect of structural dynamics, 
gravity gradient dynamics and control interactions during this effor. 
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II. DESIGN EXitMPLE 


1 . INTRODUCTION 

To understand the control problems for large space structures, a spsiiciflc 
system (the Orbiting Construction Demonstration Article (OCDA)) was selci^ted 
to develop a control system. The OCDA Is Intended to be a space base In a cir- 
cular orbit 'JO km above the Earth (Ref. 1). Its purpose Is to facilitate the 
unloading, fabrication, and assembly of objects (e.g., solar power satellites, 
microwave power transmission satellites) ferried into orbit by the shuttle. 

The OCDA has four principle parts (Fig. 2); 1) the platform or rectangular 

truss frame, 2) the boom with traveler, 3) the array that tracks the sun and 
absorbs the radiation energy using elastic membranes, and 4) the mast extending 
on both sides of the platform holding the boom and solar array on one side and 
providing a docking port for the shuttle on the other 

In its nominal attitude the OCDA is unstable with respect to aerodynamic drag 
and gravity gradient torques. These Instabilities can be overcome by active 
control that uses momentum storage devices. Environmental effects such as the 
iniip oduction of currents Into the solar array to produce magnetic control torques 
PfUi he effectively used to assist the control system in stabilising these 
;ilsturbances. We have studied stabilization by means of three reaction wheels 
with mutually orthogonal axes. 

In or 3r to calculate the control gains D In the control gain relation 
u =• Dx, the systems equations x “f(x,u) are replaced with the linear 
approximations 

X = A X + B u 

This set of equations is agumented with the cost functional 

CD 

J = u^Ru)dt 

where the weighting matrices Q, R are selected to reflect the critical nature 
of a particular node motion. The principles of optimal control then yield 
differential equations that can be integrated analytically. 

Computational expense was reduced considerably by treating the problem 
In modal coordinates for both the vibration and rigid body modes. These 











were luiiuUed in a unlforiti rather titan a hybrid winner. It was found tlwt 
satisfactory results could be obtained when only a few ef the modes were 
coiUrol led. 

In an additional study the control gains D were kept fixed while the spatial 
dimensions of the OCDA were scaled up until the point of instability was found. 


2. CONTKOi. LIMITATIONS OF CURRENT TECHNOLOGY 

As the size of a space structure increases, the structural frequencies become 
lower. Thus, for any given control system bandwidth (as determined by the dynamic 
requirements of the rigid body control specification) there exists a structure 
whoso vibration frequencies fall within the control system band. If one assumes 
that the predomlnent nuxiaj frequency of such a structure is given by the first 
mode of a simply supported beam whose dynamics are given by the Euler-Bernoulli 
equation, tlien this mode lias a frequency given by 



(In hertz) where K Is the modulu.s of elasticity, L is the length, I is the 
inertia and m is the imiss of the beam per unit length. For a slender beam the 
modal frequency is given by 



where M is the total iwiss of tho beam. The result of this analysis Is sunmuirized 
by Fig. 1, where the first modal rrequency of a slender beam is plotted as a 
function of beam lengtli. By assuming tliat this frequency is the control system 
bandwidth, we can see I'xactly bow big a structure must be for any particular 
control system specification before the rigid body control system might Infctiract 
with the structure. 

A second area where structural motion potentially impacts the control system 
is when open loop rotational c.omiiuinds are required as, for example, during a 
liigh speed slew imvneuver. The usual approach to minimizing Che bending inter- 
actions during these slew maneuvers is to tailor the shape and amplitude of the 
coiimuind pulse so that the structnrsl modes are not excited. The. optlnml control 
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Fig. 3 Structural Sice vi. Control Bandwidth. Since si7o of the structure determines 
the loweit modal frequency, the interaction the structure and rigid body 
control loop is a function of the control spec Ifi cation (bandwidth). 
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approach offers a natural framework within which such minimum structural excitation 
by command shaping may be achieved. This is shown schematically in Fig. 4 
where the command generator is shown as a linear differential equation (in this 
case the output of a filter with a maximally flat response) whose outputs are 
multiplied by "feed forward" gains. Section III-4 describes this approach in 
more detail. In either case, the exact size of a structure when slew maneuver/ 
bending coupling becomes severe enough to Impact performance depends on the control 
system pointing accuracy or figure control requirements. 

The final area where structural size can impact the control system 
performa.ce is in the coupling of disturbance induced structural motion into 
the control system. Disturbances on the spacecraft are either external — aerodynamic, 
gravity gradient, solar, thermal for example -*■ or Internal as is the case 
with vibrations induced by rotating machinery, man motions, bearing noises due to 
relative motions of spacecraft elements, etc. As an example of how these distur- 
bances may interact through the rigid body control consider the simple case 
where gravity gradient torques are used in low earth orbit to orient a long 
slender spacecraft which has a solar array mounted at one extremity (Fig. 5). 

In order to maintain the solar array at the correct orientation, the array must 
be rotated at twice the orbital rate ui . This will produce a periodic aerodynamic 

W O 

torque, due to the residual atmosphere, at the orbital period. The gravity 
gradient creates a torque which is proportionate to the rigid body rotation 
angle deviation from a line from the spacecraft t< the earth center. Thus, if 
9 is the rigid body rotational angle, the dynamic description of this inter- 
action is (for a single axis) 

18 + T^e = T sin 2w t (1) 

(jt o 

where 

Tg “ gravity gradient torque 
T = magnitude of the aerodynamic torque 
I = inertia (principle) of the axis considered 
0 )^ = frequency of an orbit 

In practice, I will be of the form I +I.sln2ta) t because of the solar array 
motion, and as such the design must 8e Basel on a time varying inertia. 
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Fig. S Disturbance Onteractions Can Cause Stability Problems as Shown by this Example 


The solution to Elq. (1) is given by 


6(t) = cos 


f 6(o) +f"“ sin /“ t 6 Co) + -t=s=- / sin / ~ (t-r) sln^ai i 

VI /t^Io Vi ° 


Now, It Is evident that If twice the orbital frequency and the gravity gradient 
frequencies are the same that 0(t) grows as t sin2to^t. Thus an unstable Inter- 
action is possible. 


For the example just described, a rather simple *'flx" is possible by 
controlling the solar array drive motor. Thus by phasing the rotation 30 that 
the solar array alternately leads or lags the sun, the aerodynamic force may be 
used to extract energy from the gravity gradient and thereby damp the gravity 
gradient control. Figure 5 shows a root locus plot that results from this 
approach which shows the gravity gradient damped regardless of the orbital 
frequency. (This design uses full state feedback for pole placement.) 
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The design uses Eq. (1) with the solar array drive provided by a motor which 
Is modeled as 


** • 

0 T 0 ■ U 

m mm 

where: 6 Is the solar array angle 

6 Is nominally 2u> 
m ■' o 

T Is the motor time constant 
m 

u Is the motor Input (u nominally is 2o) x and there is a perturbation 
du around that nominal which Is the control) 

When the solar array motor angle is substituted into the torque term on the 
right side of Eq. (1) we get 

10 + 1’ 0 = T sin0 
G m 

which, since 0 is given by 0 » u) t + A0 , can be written as 

ra mo m 

10 + T_0 “ Tsln(ii) t +A0) = T sinu) t + Tco^ tA6) 

Go o o 


Thus the dynamics of the complete open loop system (for perturbations) is 
given by 


ejsl 

U(s) 


T/I cosu) t 
o 


(S^ + T_/I)(S^ + T s') 

\J 111 


The open loop root locus for this transfer function is shovm in Fig. 5b 
plotted as the inertia I increases. The coupling described above occurs when the 
gravity gradient pole (S = ± j /T^/l) overlaps the aero pole . If a control 

system configured is shown in Fig. 5d is used, the closed loop root locus 
appears as shown in Fig. 5c) , The two zeros shown in Fig. 5c are the result of 
the feedback S + and K^S + K^. Clearly, this locus is always stable 
independent of the gravity gradient frequency. 


In general we have taken this approach for all external disturbances. 

When the optimal control model is formulated, any actuators that are available 
(the solar array drive motor for example) are included thereby allowing the 
"disturbances" to be used for active control. 

For internal disturbances, once again the magnitudes only become 
crucial if the control system bandwidth or pointing requirements are severe. 

If these sources of vibration must be damped, the structure can be included in 
the optimal control model so that they may be actively isolated. 
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To define the limitations further, part 7 of this section describes the 
analysis of the OCDA large space structure controller that we designed as the 
structure is increased in size. 

3. VIBRATION AND RIGID BODY MODES 

Before a structural analysis of the OCDA can begin the orientation and 
tilt of the solar array must be specified. "Orientation" is defined as the 
angle of rotation about the mast. The "tilt" is the rotation about the beam 
that is closest to the platform. A zero tilt means that the mast is in the 
plane of the array. Unlike the orientation, the tilt does not vary continuously. 
It is either at plus or minus 26®, depending on the longitude of the ascending 

node of the orbit. Tlie rate of orbital nodal regression is such that a flip 
is required every 22 days. The rate of radiation absorption is never degraded 
by more than 10% compared with the best tilt angle. 

The particular array position used for our study is shown in Fig. 6 together 
with the Y, Z a'fos. Their origin is at the lower end of the mast (docking port). 
The Y axis is vertical (whether up or do\m is irrelevant), the X axis is out 
of the page and parallel to the velocity. Thus the platform and array move 
edgewise. Note that the shuttle - whose structure is Ignored except for its 
moment/product of inertia - is Included. Since the shuttle is over six times 
as heavy as the OCDA, Fig. 6 indicates that the combined center of mass is close 
to origin. The two points have been treated as identical in our work. 

Each node (joint) of the OCDA is assigned a number and is tabulated with 
its coordinates and its degrees of freedom. The degrees of freedom are those 
translations and rotations about the three axes that are considered to be 
significant. Thus, up to six degrees of freedom (DOF's) can be specified for 
each node. The members (beams and strings) are also assigned numbers and they 
are tabulated with their cross sectional areas and the nodes that they join. 

Their lengths can then be found from the node table. 

If a unit force is applied at a translation DOF or a unit torque at a 
rotation DOF, the displacements for the entire DOF vector can be computed. 

The complete set obtained by varying the unit force/torque over all the DOF's 
is called the flexibility matrix. Its Inverse is called the stiffness matrix k. 
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Fig. 6 Array Oriontation and Tilt Chosen for the Structural Analysis (YZ and Y'Z' Axes 
are Shown with their Origin) 


When the dimensions (and density) of a member are given, the mass and 
moment/products of inertia can be found from simple calculations. What is 
needed, however, is the inertia matrix associated with DOF's, because that is 
how k is defined. Recall that in general there are up to six »-elevant DOF's 
at a node that joins at least two members. If only the three translation DOF's 
are significant, then the 3x3 Inertia matrix is diagonal. In the general 
case, whose theory is too complicated to be given here, for each node there is 
a 6 X 6 matrix whose elements have the dimensions of either mass, moment of 
mass, or product of Inertia. The inertia matrix for the entire structure 
will be called m. 

The original structural idealization had 1462 degrees of freedom. The 
stiffness and inertia matrices were calculated by a computer program. Computer 



11ml cations made It desirable to have a reduced order approximation which was 
accomplished using a Guyan reduction. The 249 new DOF's were a subset of the 
old except for those of Che array. There the old DOF's were translations (never 
rotations) along the XYZ axes. The new DOF's were directed along the X axis, 
normal to the array plane, and along the common normal. These directions will 
be called X, Y', Z', respectively (Fig. 6). 

New DOF's on the array were Introduced as load conditions to reduce the 
solar array model. These forces balance"should have been equal and opposite 
to the loads - this check was not satisfied at a few DOF's with the model as 
provided by Grumman structural engineering. Several attempts by us and the 
original engineer to find the source of the dlfflsulCy failed. We therefore used 
the mass and stiffness matrix as provided and Ignored the load balance require- 
ment. All this does Is to cause errors in the steady state forces on the structure. 

If X represents the displacements of Che 249 DOF's and K,M the corresponding 
stiffness and Inertia matrices, then In the absence of external forces the Guyan 
model Is: 

+ Cx + Kx - 0 (3) 

For the analysis that follows, C, the damping matrix, will be assumed to be 
zero. In Section III-6, this assumption will be corrected. 

Let us look for solutions of the form 

X “ ^ cos u)t (4) 

where u) are vector and scalar constants to be determined. After 
substituting and simplifying 

(5) 
or 

( 6 ) 

Thus 0 ) ^ Is the eigenvalue of K ^ with the corresponding eigenvector | . 

(The solution Eq. (4) Is called the vibration mode. ) 

If we Impose the initial condition k (0) = 0 there will be a 249 vector 

3 of arbitrary cons*-ants of Integration. It will be convenient to define the 

generalized coordinate vector ^ whose elements are q. cos (w, t) . The complete 

T 1/2 ^ ^ 

matrix of normalized eigenvectors j^/(^ will be denoted as $. « 


The general solution of Eq. (3) is then 

X » if 3 (7) 

If the elements of 3(0) are all zero except for q.(0) ■ 1, then x (0) 

th ^ **" 

equals the 1 vector of <I>. x(t) will remain proportional to this vector as 

its elements oscillate with angular frequency The other frequencies (modes) 

will not be excited. 

It is perhaps surprising that the matrix of eigenvectors of K~^M can 
diagonalize simultaneously both K and M. This capability is related to the 
nonorthogonalicy of the eigenvectors of which in turn is due to the nonsymmetry 

of Let and | ^ ’>e eigenvectors with distinct eigenvalues. Then from 

Eq. (5) 


0)^ M (t - K(i 
m -^m ^m 

T 

On multiplying from the left by (j)^ 


( 8 ) 


A $ » (j)^K i (9) 

m ^n -^m -n % ■ 

If these steps are repeated with the order of i , d reversed then 

“*11 

M<t> = K (|) (10) 

n -m -n -m ~n ' 

Since K and M are symmetric, the transpose of this equation Is 

A M (<) = d)’'"’ K (t (11) 

On subtracting Eq. (11) from Eq. (9) 


(o)^ - 0)^) M i =0 
m n *n ^m 


which shows that the off-diagonal elements of <I>n4> are zero and that M is 


diagonalized by 


( 12 ) 


To prove that K can also be diagonalized, the above steps 


are repeated after first dividing Eq. (8) by . 


*The eigenvectors of the matrix K”Si are not orthogonal but those of K and M 

individually are, as shown in Eqs. (11) and (12). 
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The vibration modes were augmented with the 6 rigid body translations 
and rotations that were calculated separately. The frequencies associated with 
the rigid body motion are of course zero since there are no translational or 
rotational rigid body stiffness terms in K (until the gravity gradient is added). 

The symbols K, M denote the diagonal matrices defined by 

(13) 

2 ' 

Note that Eq. (9) Implies 

The 53 vibrations modes with the highest eigenvalues (lowest frequencies) 
were computed (Fig. 7). Because the solar array consists of 11 semi-independent 
membranes, most of the modes are out-of-plane oscillations due to the first three solar 
array modes. Thirty of the latter were discarded by retaining only one of these 
membrane modes per frequency leaving 23 vibration modes. The corresponding 
terms of K, M were also deleted which corresponds to reducing the order of the 
system by assuming the frequencies are zero. The exact mathematical statement 
of this "singular perturbation" is the followin|: 

If the modal coordinate vector is denoted by g then g is partitioned into 
components g^^ and g 2 where. In terms of this partition 


^^11 

0 


ft ^ 

3l 


1 

(-• 

M 

0 ' 


r 1 
3l 

0 

M 22 




0 

^22 

m* 


32 

L. mJ 


where all of the matrices are diagonal (as in Eq. (13)). Now If M 22 has a 
constant e factored out of each of its terms, then (see Ref. 2), the "reduced 
order" model becomes g^ + g^^ = 0 as £ -*• 0 (i.e., g 2 0) . This tacitly 

implies that the full 249 vector x is given by (7) with g 2 ® 0. 

That is 


ss 

~<t> <I> ” 

11 12 


'3i' 

B 

*n 3i 


A d) 

[21 ^22 


0 


.*21 31. 
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FREO 

He 

LB, 

SEC^/IN, 

DESCRIPTION 

•1 

0.04766 

10,653 

SOLAR ARRAY - 1ST NORMAL TRANSLATION 

•2 

0 05538 

3.913 

SOLAR ARRAY • 1ST TORSION (1ST MEMBRANE) 

•3 

0.06008 

3,822 

SOLAR ARRAY - PANEL MODES (1ST MEMBRANE) 

4 

5 

0,06135 

0.06176 

3,813 

3.475 

SOLAR ARRAY - PANEL MODES (1ST MEMBRANE) 
SOLAR ARRAY - PANEL MODESOST MEMBRANE) 

6 

0.08192 

3802 

SOLAR ARRAY - PANEL MODES (1ST MEMBRANE) 

7 

0.08201 

3 711 

SOLAR ARRAY ~ PANEL MODES (1ST MEMBRANE) 

8 

0.06206 

3,743 

SOLAR ARRAY - PANEL MODES (1ST MEMBRANE) 

9 

0,06209 

3.647 

SOLAR ARRAY - PANEL MODES (1ST MEMBRANE) 

10 

0,08211 

3,787 

SOLAR ARRAY - PANEL MODES (1ST MEMBRANE) 

11 

0.06212 

3,697 

SOLAR ARRAY - PANEL MODES (1ST MEMBRANE) 

12 

0.08213 

3,736 

SOLAR ARRAY - PANEL MODES (1ST MEMBRANE) 

13 

0.06214 

3,676 

SOLAR ARRAY - PANEL MODES (1ST MEMBRANE) 

•14 

0.08400 

1 1 .240 

SOLAR ARRAY ~ 2ND NORMAL TRANSLATION 

•16 

0.09226 

20,274 

SOLAR ARRAY - SIDE BENDING (X) 

•16 

0.10329 

6.128 

SOLAR ARRAY - TORSION (2ND MEMBRANE) 

•17 

0,10961 

28,561 

HEAT TORSION - BOOM + PLATFORM AND ARRAY LOCATION 

•18 

0.11288 

3.625 

SOLAR ARRAY - PANEL MODES (2ND MEMBRANE) 

19 

0.11404 

3,754 

SOLAR ARRAY - PANEL MODES (2ND MEMBRANE) 

20 

0.11442 

3.694 

SOLAR ARRAY - PANEL MODES (2ND MEMBRANE) 

21 

0.1 1459 

3.767 

SOLAR ARRAY - PANEL MODES (2ND MEMBRANE) 

22 

0.11468 

3.699 

SOLAR ARRAY - PANEL MODES (2ND MEMBRANE) 

23 

0.11473 

3,754 

SOLAR ARRAY - PANEL MODES (2ND MEMBRANE) 

24 

0,11477 

3,704 

SOLAR ARRAY - PANEL MODES (2ND MEMBRANE) 

25 

0,11479 

3,753 

SOLAR ARRAY - PANEL MODES (2ND MEMBRANE) 

26 

0.11480 

3.682 

SOLAR ARRAY - PANEL MODES (2ND MEMBRANE) 

27 

0,11481 

3.746 

SOLAR ARRAY - PANEL 10DES (2ND MEMBRANE) 

28 

0,11482 

3,715 

SOLAR ARRAY - PANEL MODES (2ND MEMBRANE) 

•29 

0.11628 

26,743 

riTICAL BOOM + PLATFORM SCISSOR MODE 

•30 

0.12259 

8,229 

; 0,,AR ARRAY - PANEL MODES (3RD MEMBRANE) 

•31 

0.13057 

5.362 

... i-AR ARRAY - PANEL MODES (3RD MEMBRANE) 

•32 

0,14846 

3.223 

SOLAR ARRAY - PANEL MODES (3RD MEMBRANE) 

33 

0.14961 

3.926 

SOLAR ARRAY ~ PANEL MODES (3RD MEMBRANE) 

34 

0,14984 

3.425 

SOLAR ARRAY ~ )’ANEL MODES (3RD MEMBRANE) 

35 

0.14992 

3,791 

SOLAR ARRAY - PANEL MODES (3RD MEMBR.9NE) 

36 

0.14996 

3,747 

SOLAR ARRAY - PANEL MODES (3RD MEMBRANE) 

37 

0,14998 

3,744 

SOLAR ARRAY -- PANEL MODES (3RD MEMBRANE) 

38 

0.14999 

3 805 

SOLAR ARRAY - PANEL MODES (3RD MEMBRANE) 

39 

0,15000 

3.777 

SOLAR ARRAY - PANEL MODES (3RD MEMBRANE) 

40 

O.15C01 

3,724 

SOLAR ARRAY - PANEL MODES (3RD MEMBRANE) 

41 

0.15002 

3,683 

SOLAR ARRAY - PANEL MOucS {gFO MEMBRANE) 

42 

0.15002 

3.6033 

SOLAR ARRAY - PANEL MODES (3RD MEMBRANE) 

•43 

0,15178 

7,6922 

SOLAR ARRAY - PANEL MODE 

•44 

0.15395 

4,6886 

SOLAR ARRAY - PANEL MODE 

•45 

0,17424 

16,5774 

PLATFORM NORMAL BENDING 

•46 

0.20019 

3.24636 

SOLAR ARRAY - PANEL MODE 

•47 

0.22621 

48.6614 

PLATFORM TORSION 

•48 

0.25980 

5.6988 

SOLAR ARRAY - PANEL MODE 

•49 

0.26885 

10,4788 

SOLAR ARRAY - IN-PLANE 

•50 

0,27068 

7,5181 

SOLAR ARRAY - PANEL MODE 

•51 

0,30673 

4,1957 

SOLAR ARRAY PANEL MODE 

•52 

0.30795 

210,487 

PLATFORM LATERAL BENDING 

•53 

0.32406 

200.3 

PLATFORM BENDING 


•THESE MODES HAVE BEEN RETAINED 

0357 007W 
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4. EXTERNAL FORCES AND TORQUES 


If an external force/torque vector f is present » Eq. (3) generalizes to 

M X K X - f (14) 

Equation (7) can be used to reduce the number of equations from 249 to 29 by 
substituting x " in (14) 

M <l)g K - f (15) 

T 

Multiplying on the left by ‘I’ and using Eq. (13) 

Mg Kg - *'^f (16) 

This equation is exact although x as calculated from Eq. (7) is approximate. 

The complete orthogonallzatlon procedure described above follows the 
classical structural analysis techniques. A more general approach utilizes the 
linear algebra of symmetric matrices. Thus (14) can be rewritten as follows: 

• Let M be written as a product of a lower and upper triangular factor as 

M - LL^ (17) 

This "Cholesky" factorization may always be performed since M will 
always be symmetric and positive definite. 

• Define a new vector z as 

z ‘ X (18) 

• With z replacing x In Eq. (14) we get 

z = -l’^KL"^z f L“^f (19) 

-1 -T 

• Modal Transformation: Since K was symmetric L KL remains symmetric 

and the following results from linear algebra may be used to diagonalize 
It. 

a Lemma (Ref. 3): A real symmetric matrix can always be diagonalized by 

an orthogonal transformation and its eigenvalues are always real. 

T T 

Thus, let z = <I>'g' where il>* Is orthornormal (<I''iI' ' <!► «>* » I) so that 
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where 



and we assume < 


4 ^ 




Notice that the 51' defined In Eq. (20) and the g from Eq. (16) are not Identical 
since M In Eq. (16) Is not the Identity matrix. In fact 3' » M g and 
$ a MO' where 



The only external generalized forces considered were the aerodynamic drag 
acting on the solar array and the external torques due to the gravity gradient 
and the moment \im wheels. 


Before discussing the drag we digress to recall that the axes used for the 
earlier computer work that was taken over for the present Investigation were those 
of Fig. 6. However » the earlier report (Ref. 1) followed the more conventional 
usage and Interchanged the labels for the Y and Z axes. When the solar array 
has the orientation shown in Fig. 6 it Is moving edgewise and there is no aero- 
dynamic drag. We used the axes of Ref. 1, but in order to Investigate 
stabilization in the presence of drag we rotated the OCDA so that the array 
moved nearly facewise (Fig. 8). Note that motion Is along the Z axis and the X 
axis Is vertical. 
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since Che drag was distributed over the solar array membranes* It was 
necessary to find lumped values at the nodes. The drag at a corner was assumed 
to be one-fourth that of an Interior node. The drag at other edge points was 
one-half that at interior nodes. The forces were assumed to be In Che direction 
normal to the array (-Z*, Fig. 8) rather than opposed to the velocity as they 
should have been. It was impossible to give the drag the latter direction because 
Che Y* and Z* DOF's were not selected In pairs at each membrane node. Indeed 
the latter were always selected but never the former. 

These remarks will become clearer if we look at <t> in more detail. It has 
the form 


m 



12 




0 



0 



( 21 ) 


where in this representation the first column stands for the 23 vibration modes, 
the second the three rigid body translation modes, and the third the three 
rotation modes. The 3x3 unit and null matrices are denoted by I and 0, 
respectively. The matrices M, K and the vectors f and g can be decomposed 
similarly into 


“ dlag ... M22,23^» ^2 “ ^”24, 24 ^*26,26^ 


M3 = dlag tM27^27***^29,29^» ^1 “ ^^1,1 •** ^23,23^ 


( 22 ) 


K2 = K3 = 0 


.T T r T . T, 

- “ * £2 ’ ^3 ^ 


T T T T 

3 “ ^3i > 32 • 33 ^ » ^1 ^2 33 


The nonzero elements of the 243-vector f^ are the drag terms discussed above, 
the three-vector I2 is null, The three-vector £3 of gravity gradient and 
momentum wheel torques will ba discussed in detail below. Equation (16) can 
now be written as the three equations 
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( 23 ) 


3l + .002(Aiij)‘‘ ij + Ki 3l - hi h * hi ^3 

”2 32 " *12 ’ ”3 33 " *13 I). * 

Note that a modal structural damping term has been introduced into the vibration 
modes. This term changes the solutions of the homogeneous equation (Eq. (3)) by 
a factor exp(- .OOloit) . It also has a major effect on the modes Which were 
reduced out of the dynamics as will be discussed in Section 111-6. 

The second of Eqs. (23) represents the three equations 

”2^ “ ^24-^1 ’ ^2^ * ^5^1 ’ ^2^ “ ^26^1 

where ^2 * M 2 ^ 24 “ ^25 25 * ^26 26' ^ coordinates of the center 

mass; and (1 =» 24, 25, 26) is a 243-row vector whose elements equal the first 

^ th 1 

243 elements of the 1 vector of 4>. The terms of the scalar product 

are all zero because the nonzero terms of ^24 never coincide. The 

terK(s of ^25 that coincide are all equal to (using the axes of Fig. 3) 

cos (Y,Z’) = sin 26°, the corresponding terms of are cos (Z,Z') = cos 26°. 

The acceleration In the Y direction does not vanish although It Is normal to the 

velocity. There is no way to correct this qualitative error because no Y* DOF's 

were selected for the array membranes. 

Reference 1 states that the momentum wheels are located in the mast midway 

between the platform and the boom. However, they were placed at the docking 

port (origin) for the present study. We have retained the latter position and 

assumed that the angular momenta h ,h , h of the three wheels have fixed 

® wx coy u)z 

directions along the coordinate axes. Small control torques u = h , 

^ X cox 

u = h , u = fi are obtained by changing the magnitudes of the angular 

y coy z coz ^000 c 

momenta. In accordance with conservation laws, an Increase in the angular 
momentum of a wheel produces a decrease in the angular momentum of the OCDA. 

The rigid body yaw about the X axis (Fig. 6 ), pitch about the Y axis, and 
roll about the Z axis will be denoted by a = q^y, p = y = q 2 ^, respectively. 

The corresponding moments of inertia are A - M 27 27 * ® ^28 28’ ^ ~ ^29 29* 

Linearized gravity gradient torques are given In Ref. 4, p. 244 and appendix 1. 
Thus the detailed expressions for the last of Eqs. (23) are 
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(25) 


A“ - ^27il + 

»6 - Ijsil + -Wz - “y 

Cy ■ "liirtf, + ay +bA “W h + w h - u 
' *29-1 y wx X wy z 

These equations use the following abbreviations : u)^ is the angular frequency 

corresponding to the 90 minute orbital period » 

a = 4w^(A-B), b = o)^(A-B+C) 

c » a)^(C-B) , e - 3oj^(A-C) (26) 

o o 

tt)*d + wy, a)=^-u) ,u)"y - o) a . 

X o y o z o 

The scalar products ^28-1’ ^29-1 zero. This occurs because 

DOF's were selected in a synunetric manner over the array so that there is a 
negative term to cancel every positive term. Thus defining the drag to be in 
the Z' rather than the -Z direction does not introduce qualitative errors in the 

torques as it does in the forces on the center of mass. The quantitative 

error for the torque about the X axis tends to cancel out when summing over 
all the DOF's of the array. That is, the moment arm is smaller for Z' forces 
than it is for -Z forces at nodes near the platform and larger at nodes far 
from the platform. Note that another non-homogeneous term enters through 

Equations (23) with initial conditions £(0) = g(0) = 0 describe a step 
disturbance - the drag is zero when t < 0 and jumps to its full value when t “ 0. 

It is also interesting to study an impulsive disturbance occurring at t “ 0. 

In this case the non-homogeneous terms are removed from Eq. (23) and the initial 
conditions g(0) are set equal to them. 
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5. STATE VARIABLE MODEL FOR (X)NTROL DESIGN 

The optimal control model we arc developing uses 23 modes by leaving out the 
less Important vibration modes through order reduction In modal coordinates. 

The force vector f consists of a distributed aerodynamic load on the solar 
array and the gravity gradient torques on the rigid body (we assume that the 
gravity forces are distributed as described In Appendix I). Thus the force vector 
f becomes 


h 




243 vector of aerodynamic forces distributed 
Into the structure 


3 vector of gravity gradient torques and 
control torques 


(27) 


unere 


-3 


Hu 


T- + T^ X. + H u + H as In Eq. (25) 
U 2 —j ■“ o 


control forces: H “ 

’ o 


’tc h 

O 0)Z 


0) h 

O U)X 


is the non-homogeneous forcing function. 


T = gravity gradient torques Induced by rigid body rotations 
^1 

Tp = gravity gradient torques Induced by rigid body rotational rates. 


Following Eq. (25) T and T^ are given by 


% (C - B) 


3o)^ (A - C) 


4o)^ (A “ B) 
o 


(28) 


23 




Hu 

mm 


0 

0 



u»^(A-B+C) 

0 

0 


0 

> 

1 

0 


0 

N 

3 

1 

O 

0) 

y 

0 



“>2 0 

-0) 

X 

0 


^wy 

-w w 

y X 

0 

ilD 


Kz 





_ Array, 


( 29 ) 


(30) 


h - angular tnomentum of the controller about the principle axes 

^Xf V f 2 


)ID6 Array - torque produced by the solar array angle and £ is the 

lever arm of the drag force D. When D Is non-zero f^^ is 
zero and vice versa 


As can be seen the formulation of ':he control terms Eq. (30) allows the 
simultaneous incorporation of the gravity gradient, aerodynamic drag and 
control momentum exchange devices. To distribute the terms from Eq. (25) 
Into the structural mode model Eq. (16) we let f be given by Eq. (27) In 
Eq. (16) and get 




~ 

r 

1 . ~ 

Jl 

0 

0 

J 

3i 

0 

^2 

0 


32 

0 

L 

0 

M 

3j 

1 

i3_ 


0 0 


’ 31 “ 

N 

0 0 0 


32 


o 

o 

o 

1 


33 

> 


r.T ^ 


r ^ 

^31 



.T 

^32 

< 

’^G/*31*32»33i a + \ '*31»32*33> ^ \ 

*T 



®33 


< J 


( 31 ) 


2A 


which niay be reduced to give 


S 


.-1 




- q : + M M 


4>- 


31‘G/31 1 31*G, 32 "1 '^31*G, 33 


^2^^32’’g^^ 31 ^2 "^32Tg^^32 ^2*^32'^G^'®*33 


,-l.T 


-l.T 


-IT -IT 


-1 T 


‘3 33‘G/31 "3 "33*G, 32 “3 "33^G, 33 


+ 


”.-l T „ -1 T „ 1 


r. -l.T 

M T il) ... M 4 T 

1 *31 Gj^31 ••• *^31^G^*33 


M 4> 

1 31 

..-l.T „ . w-l.T „ . 


.-1 T 

^2 '*’32^0^*31 •" ^2 “^32^Gj^^33 

2 -^ 

M iJ> 

“2 *32 

..-1*T . -l^T 


-1 T 

M $ T if . . M 7 ^ 

3 33^g/31 ••• ”3 33 G, 33 


^<^33 

^ 1 1 J 


J 


V / 


Hu+M”^H 

o 


T 


(32) 


Y <p • ^ 

and finally, by defining the state vector z as z " (q : q ) the model 
used for designing the optimal control system becomes: 


r- 

0 

0 

0 

cn 

CM 

:3 

0 

0 


0 


“0 

0 

0 

0 

0 

^3 

0 


0 


0 

0 

0 

0 

0 

0 

"3 

2 + 

0 

— I T 

u + 

0 

^11 

®12 

^13 

T 

11 

T 

12 

"^13 




M H 

*2. 

S21 

S22 

S23 

'*^21 

T 

22 

"^23 


-1 T 
''2 ^ 32 « 



S 31 

S32 

^33 

^31 

T 

^32 

T 33 


-1 T 














or neglecting the M term 

z = Az + Bu 


(33) 


(34) 
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where 


2 is now a 58 vector (46 states for the structural dynamics and 12 rigid 

body states) 

S and T are defined Implicitly by Eq. (33) 

Is the n X n Identity i^natrlx. 

The nominal solar array angle is 23“ so that the last component of u acts as a 

-IT 

disturbance on the structure through the ^ 31 ** component of the matrix B In 
Eq. (34). Thus, the uncontrolled motion of the spacecraft can be developed 
and used as a reference to compare the controlled motion. The open loop dynamic 
response at various times are shown in Fig. 9. For the control system design 
in the next section, the rigid body tranlatlon modes have been deleted so the 
control vector z Is reduced In order to 52. 

Notice in the formulation of the dynamic model leading to Eq. (33) that 
the rigid body modes are retained in the same coordinates as the structural 
modes. This is distinct from Likens et al. (Ref. 5) where a hybrid coordinate 
system is used. 

6. OPTIMAL DIGITAL CONTROL FOR OCDA 

It is not the main point of this report to present the theory of linear 
discrete-time optimal control. Therefore, t^e will briefly outline here the 
major assumptions made in designing optimal discrete systems. For further 
details the reader is referred to the excellent tutorial paper by Dorato and 
Levis, Ref. 6. 

The continuous model of the spacecraft in state variable form (where all 
matrices are constant) is given by Eq. (34) 

^ = Ax + Bu + Cw 

jr = Kta + V (34) 


The usual notation in the control literature is to use x as the state vector. 
This convention is used here thus the original version of Eq. (34) (on page 25) 
is modified (in the original formulation x was the physical motions of the 
finite element degrees of freedom) . 
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Fig. 9 UncontroOled, Open Loop, Response of the OCDA — impulse Loading of the Soler 
\ ~ Array as the Driving Force 

\\ 
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Fig. 11 Orbiting Conitant Biia - No Controt - Modai 6, 6 and 7 Unchn) 
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Fig. 12 Orbiting Constent Base — Open Loop Rigid Body Rotation Modes (Radians) 
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where x ie the 52 vector of spacecraft states and any augmented noise states 
(which are not Included) u Is a 4 vector c^outrol; w would be a vector noise 

process (which Is assumed to be a white process); )r Is the vector measurement 
which we are assuming Is x; and v Is an m 2 vector white measurement noise. 

The discrete-time version of Eq. (34) Is given by 

Xjj^^ - 'I' (At) Xj^ + r(At) + G(At) Wjj 


*K 




At " the sample time 

AAt 

<^(At) " the transition matrix evaluated at At (l.e.» e ) 

At 

r(At) - / ' 4(At-T)B dr 
o 

(i.e., f - Ar+B, r(0) - 0 r(At) » / e^^Bdt) 

o 

T 

G(At) * a matrix which satisfies G(At)G (At) “ P(At) where 

P(At) Is the covariance of matrix of x from Eq. (34) at t » At. 

In general, one would like to solve the following problem: "Given a 

quadratic performance Index of the form 

T 

j - f (x^Qx + 

t 

o 

find the control which minimizes J subject to the constraint of the 
differential Eq. (34)." 

The discrete version of Eq. (35) is obtained from Eq. (35) and the fact 
that the controls are held fixed over the sample Interval At. This discretization 
of Eq. (35) Is necessary to be able to compare the discrete design with the 
continuous design. It Is generally better to use the performance Index for 
the discrete system which matches the continuous system also to determine the 
effect of Increasing sample time. Thus, the performance measure Eq. (35) becomes 

J = Zx^Qx^ + 2uJsXj^+i^Rt^ (36) 


31 


where 


Q - / ^^(t) Q ^(t) dT 

o 

At 

S - / r^T) Q <Kt) dT * 

o 

At „ 

R - / {R + r^(T) qr(T)} dr 

o 

The computation of Q, S and R uses an eigenvalue-eigenvector approach where the 

At -1 

transition matrix <I>(t) Is written as Te T with A as the diagonal matrix 
obtained by dlagonalixlng the matrix A. 

The discrete-time control problem Is solved by using either the discrete 
maximum principle or dynamic programming. In either case the ultimate solution Is 
that ^ Is a linear function of the states of the system where the gain mtrlx 
satisfies 

- - (R + + s) (37) 

and solution of the discrete matrix Rlcattl equation 

- + s)'‘'(R + + s) (38) 

with Pjj = [0]. 

If, as Is generally assumed, the steady-state (constant) gain matrix is 
desired, we have found that the best way to solve for the gain, Eq. (37), is 
to use Potter's technique (Ref. 7) of evaluating the eigenvalues and eigen- 
vectors of the 2n x 2n matrix which results from the application oi the discrete 
maximum principle. Problems such as the OCDA control which are high order and 
with significant differences in time constants are as easily solved as low-6rder 
pr:-.,lems when using a numerical eigenvalue, eigenvector technique. If such 

The matrix T(t) Is defined after Eq. (35); also It Is^signif leant that the 
discrete performance Index contains a nonzero matrix S even though Eq. (35) 
does not. 
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problems were solved using Iteration, the computer time required would be 
exorbitant. This matrix, for the performance index Eq. (36) is given by 


■ “11 

« 12 ‘ 

_ «J 1 

«22 

w 


where 


”ll " “ rR‘^S)"^ 

“ (4* - rR“^s)"^rR“^r^ 

- (Q S^R S)(* - rR“^S)"^ 


T ~T~-1 T ~ -T~-1 "* 

H22 - (® - s^R ) + (q - s^R ^s) 


- I’R ^s) ”^rR r 


(39) 


H has the property that its eigenvalues are the n poles of the closed-loop system 
and. their reciprocals (there are 2n eigenvalues, n stable and n unstable). 

The steady-state gain is given by the Potter algorithm as follows; 

Let W be the nmtrlx of eigenvectors, then 



where A represents the eigenvalues outside the unit circle. Then if W Is 
partitioned as 


■ “n 

"12 ] 

"21 

5s 




the steady-state cost Is given by 



is n X n 


P 

o 


”21^11 


-1 


(40) 
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The program that Is used to derive the optimal controls uses this technique 
to obtain the optimal gains. The Q - R algorithm Is used to determine the eigen*- 
vector matrix W. Potter's technique applies to both continuous and discrete 
systems; however, in the continuous version the eigenvalues of the 2n x 2n 
matrix H have symmetry about the imagiaary axis (rather than radial symmetry) . 
Since we use the continuous optimal design as a reference, the Potter solution 
for the continuous design existed first. To take advantage of this existing 
program we map the roots of the matrix H using the bilinear transformation 

w * (1 - z) / (1 + z) 

This allow.s the roots Inside and outside the unit circle to selected based 
on their locations in the left-half or right-half w-plane. 

The last point we have to consider is how one selects the sample time At. 

This has already been discussed in a former paper by the authors (Ref. 8), but 
a brief outline Is in orc’er here. The selection of a fundamental sample time 
uses the fact that between samples the system is essentially open- loop. 

Therefore, If perfect control at the sample times were achieved, between 
samples the uncertainty would propagate via the covariance matrix differential 
equation 

, m m 

P = AP + PA + CC 

P(0) = [0]. (41) 

At the first time t that the uncertainty due to this open-^loop propagation 
exceeds a bound specified by the control specifications, there must be a sample 
to lower this uncertainty. Since we are neglecting uncertainties that are due 
to the control (feedback) and the state estimation in doing the propagation of 
Eq. (41) (l.e., P(0) is really not fO]), the actual sample rate should be 
modified based on the closed-loop noise response. As we will see the existence 
of computer aided design programs allows all of the above steps to be implemented. 
To show all of the output of the computer computations for the OCDA example 
would not be very instructive; we therefore have shown In L^ection IV an exercise 
for a simplified example. Reference 9 shows more details of this analysis. 

Many times It is necessary to provide a constant force to provide the desired 
steady state position. To guarantee that this force does not cause an undesirable 
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deviation at some position on the structure, the control is assumed to be applied 
through a set of Integrators. Thus f “ Bu’, where B dlstrubutes the control 
forces u' through the structure, u* is integrated control, l.e., 

t 

u’ » / u dt 

o 

or 



The dynamic state of Eq. (34) is modified as follows: 
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(42) 


is the ’’state variable" model of the structure which is used to design the control 
system. ^ has dimension 52 + p where p is the dimension of u (the control) . 

The optimal control design based on the model Eq. (34) uses the performance 
measure Eq. (35). But Eq. (35) is really based on the actual full 243 degrees 
of freedom and through the order reduction, Eq. (35) becomes 
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Equation (43) must be used Instead of Eq. (35) to develop the discrete control 
performance measure Eq. (36). 

As can be, seen a cross product between and u is now explicitly Included 
in the performance measure. Furthermore, even though is a reduced state, 
its order is still quite large. The solution of Eq. (43) as an optimal control 
problem requires extrismely accurate nianerical techniques. Furthermore when 
the solution has been achieved, there is still the problem of the control exciting 
the higher frequency modes. Notice that this design will require full state 
feedback (l.e., 3^) • This implies using measurements of m physical degrees 

of freedom, x*" and x*" . Then since 

g = <?L X 

(44) 

3l = '■'ml S” + V 


36 


T 

where L denotes the restriction of the matrix $ L to an mxm matrix corresponding 
m 

to those nodes x"* measured. No Kalman filter or observer is Introduced by 
implementing this pseudo-full state measurement scheme. 

The determination of requires some explanation. In the order reduction 
described in Section II-3, the fact that the "fast" states (denoted by ^ 2 ) 
were zero meant that . 


r a 

*11 3i 

- '**21 3l_ 


(45) 


T 

where was the m vector of retained modes. Since 4> ^ “I (.as formulated in 
Eq. (20)), the full mode vector g is given by 
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where x^^ is an arbitrary subdivision of the node degrees of freedom whose 

dimension is the same as g.. Then, since our assumption for the order 

„ 2 -1 T T 

reduction was that g 2 - which implies g 2 = (5^2^ ^‘^12'^22^ I “ 


*12 2^1 *22 2f2 ' 


(47) 
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therefore 


?2 ° ^ *12 22l " “£2 
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Thus, if is the measured degrees of freedom x , and are given by 


•'n, ' 
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The optimal control design for the structure defined above was developed In 
a sec^uence of four different designs* as follows: 

• Rigid body control only; the resulting controller is then 
tested using the complete structural dynamic model to 
verify the rigid body control stability 

• A continuous optimal control using only the momentum storage 
devices for control 

p A discrete optimal control as In case 2 

• A control design as in cases 2 and 3 which uses only a reduced order 
model; the resulting controller Is again verified using the full 
dynatrdc model. 

The table below summarizes the results for each of the designs. 

RESULTS OF OPTIMAL CONTROL DESIGNS 
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As can be seen the gravity gradient is unstable in roll. To achieve the level 

-3 

of damping shown (critically damped rigid body and from 2 x 10 to 0.4 on 
various flexible modes) the control forces were not excessively large. For a 
typical response (the same as the open loop conditions shown in Figs. 9, 10, 11, 
and 12, where an Impulsive force is applied to the solar array) the control 
torques are as sho^ in Fig. 13 (continuous control) and Fig. 16 (discrete control) . 
As can be seen, no more than 8 ft. lbs of torque are required to control the 
vehicle and to damp the flexible mode vibrations created by the 20 ft. lb. 
aeroydnamlc torque. The time response using both the continuous and discrete 
control at .2 second sample time are shown in Figs. 14 through IS and Figs. 

17 through 18. The time histories of the complete structure are shown in Figs. 

19 and 20 for continuous (Fig. 19) and discrete control at 1 second sampling 
time (Fig. 20) . 

The assiunption that full state feedback is available is not overly restric- 
tive. It is quite easy to Implement this concept using strain gauges or an 
optical system for sensing the motion of the structural nodes. Since, for 
most applications, the noise equivalent strain is so low filtering is not 
necessary which eliminates the sensitivity of the control to the’ modal parameters 
Introduced by a Kalman filter or an observer. A Kalman flltei; by Introducing 
notches at the modal frequencies, causes a sensitivity to the mode frequencies 
that outweighs its beneficial effects for the spillover of measurements from 
the truncated modes (see Section III-3). 

7. EFFECT OF STRUCTURAL SIZE CHANGES ON CONTROL 

We have investigated the stability of the OCDA space structure control 
system as its dimensions are rescaled while the feedback ''control matrix is held 
fixed, which corresponds to a spacecraft i\7ith Increasing dimension using a fixed 
gain controller. The feedback matrix D is first calculated so that the system 
X = Ax + Bu, (u = Dx) is optimal where x is the 52 vector of modal displacements, 
rates and rigid body angles, and rates. Then A, B are replaced by matrices 
corresponding to a structure all of whose linear dimensions are scaled by a 
factor L. We seek the value of L at which the system becomes unstable. 


These figures are frames from a 16mm sound movie that was produced during 
this contract. Copies of this film are available. Interested parties 
should contact the authors. 
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Fig. 16 Orbiting Conitruetion Bate - CloMd Loop Oiicrete (0.2 Soc Simpit Tima) - Control Torquai (Inch Lbs.) 
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Fig. 18 Orbiting Construction Bate - Sample Time 0.20 — Rigid Body Rotational Modes (Radiant! 
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Fig. 19 Controlled, Continuous, Response of the OCDA 
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Recall that the OCDA is gravity gradient unstable and that D is determined 
by the rotations at the center of mass. The aerodynamic disturbance was assumed 
to be an impulse rather than a step. Our analysis employs 23 flexible modes. 

There are 18 translational modes whose generalized masses are scaled by a factor 

3 

of L as though they were true masses. The five remaining modes are rotations 
and their generalized masses correspond to torsion and are scaled by a factor of 

5 

L , as though they were moments of Inertia. This is also true of the rigid body 
rotation modes. The elements of the stiffness matrix are similarly multiplied 

3 

by either L (bending) or L (torsion). The three vectors of the mode matrix <1’ 

that correspond to rigid body translations do not change because their nonzero 

elements are dimensionless. (We have been deleting these from our analyses 

in any case.) The linear nodes of the flexibility and rigid body rotation modes 

are scaled by L while the angular nodes are kept fixed. Note that although the 

2 

impulsive aerodynamic drag Increases by L , it plays no role in this analysis 
since it affects only the initial conditions for the state x. The gravity 
gradient torque is also affected by the scaling, but since these terms appear 

3 

in the stiffness matrix they are automatically scaled by L . 

Some root loci are shown in Fig. 21. It can be seen chat instability 
occurs at L = 2.4 (Fig. 21a) and for the design which Ignores the flexible motion 
(rigid body only), the system is stable for L up to 5.2 (Pig. 21b). In both 
cases the instability is caused first by the rigid body/gravlty gradient going 
unstable. The structural mode that goes unstable is the same in each case but 
occurs at a larger L for the case of full state feedback. Thus we can conclude 
that the optimal designs are reasonably insensitive to variations in the structural 
size. 

8. RIGID BODY CONTROL FOR PRECISION POINTING 

The problem of precise control of the vibration of a system and the control 
of the rigid body are incompatible problems. To control vibration requires 
a wide band system while the usually high noise on the attitude sensor relative 
to the rate sensor tends to make the bandwidth of the attitude loop very small. 

This suggests that a hierarchy of control loops should be designed - the first 
loop a narrow band attitude loop, the second loop a wider band rate loop that 
controls the rigid body rate and the final loop a rapid moving (wide band) loop 
that achieves vibration control. We describe in this section how the first of 
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Fig. 21(b) Optimal Rigid Body Control Design - Root Locus As Size (U Increases 
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these loops - the attitude control loop can be designed around the second rate 
control loop. First, however, we show that the rate sensor noise relative to the 
attitude sensor noise determines the bandwidth of the attitude loop. 


\Jhen sensor noise and external system noise are both exciting the control 
loop, the design of the control system requires that an estimator (a Kalman 
fil<*er) be designed to provide the best estimate of the rotational rigid body 
states. Consider the rigid body dynamics given by 
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X 


+ 


[ 




w 


(50) 


jr “ 5 S 


that is, the rigid body dynamics are 10 » noise where in Eq. (50) Xj^ “ 0, 

X 2 ■' 0, 1 is the inertia of the rigid body axis and the measurement ^ is a 
measurement of attitude and rate with additive noise n. We assume that the 

2 2 ”2 

noises w, nj^ and n 2 are white with variances and 02 , respectively. 

The continuous time Kalman filter is given by 

i fo ll . 

X = X + K(x - x) (51) 

[o oj - 

where K « P 
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The significant aspect of the filter is what it does to the control system 
when it is inserted in the control loop. It is well known (See Ref, 7) that the 
closed loop system with a filter in the loop has closed loop poles idiich. are those 
of the optiiiiial system cascaded with those of the filter. Thus, the bandwidth of 
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the closed loop system will be that of the optimal control if the filter band- 
width is greater than the control. If the filter bandwidth is smaller than the 
optimal control, then the closed loop response is dominated by the filter. To 
determine the poles of the filter (Eq. (51)), we use Potters method on Eq. (SO). 
Thus, the filter poles are the left hand plane eigenvalues of the 4x4 matrix 

0 0 +l/o^ 

+10 0 

0 0 0 

0 o^/I 0 

w 

These eigenvalues may be determined analytically as the solutions of the charac- 
teristic polynomial of Eq. (52). This polynomial is 
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where and A 2 are the left half plane poles (with the minus sign). To 
characterize the bandwidth of the filter, the following parameters are introduced. 
Let (the signal to noise ratio for the rate sensor) be 0 Jl/o 2 and let p 
(the ratio of the attitude sensor noise to the rate sensor noise) be 


Then 


A^ = - /2/2 { 


A2 = - *^/2 { 


1 + /iT^ 


4/(S^p) } 


x./TT 


1/2 


1/2 


(54) 
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Figure 22 shows the root locations for various conditions of and p. The 
situation in Fig. 22a is the critical one. It describee the case when p is large, 
i.e., when the attitude noise is larger than the rate sensor nols. For a typical 
precision rate gyro and star tracker, the ratio p is on the order ri 100 which 
is clearly the case of Fig. 22a. Thus, the bandwidth for the attitude control 
loop is determined by the filter pole near the origin. In other words, the 
settling time for the attitude loop is determined by the time it takes to 
determine the attitude from the rate measurement. Since this is the case, 
an attitude control system should be designed that uses rate sensors for the 
basic control and the attitude control loop should provide rate commands. When 
this is done it becomes possible to also estimate some of the rate gyro parameters. 

The basic function of any onboard attitude control system is to maintain 
constant LOS of the spacecraft axes. Ideally, it should be free of long-term drifts 
and jitter, and steerable by ground command. To determine the degree to which 
these ideal characteristics can be achieved, and by what techniques, a baseline 
approach to the control, determination, and parameter calibration problem is 
proposed and analyzed. By considering a generic design such as the one proposed, 
it is possible to identify the specific sources of LOS error and to propose and 
evaluate techniques for their compensation. 

The proposed design utilizes rate stabilization with a three-axis strapdown 
gyro package coupled with a set of coarse and fine reaction wheels. The 
stabilization loop recieves ground commands of the required rates it must 
hold constant, with respect to inertial space, and the period over which they 
must be held. This system is an onboard velocity feedback loop which can 
respond rapidly and the position loop is closed through the ground, using 
an onboard star sensor. The position loop can be closed through a ground loop 
because of its very low bandwidth as described above. Figure 23 depicts the 
essential elements of the proposed approach. The dashed portion represents the 
onboard stabilization loop, the remaining portion represents the ground-based 
position loop. 

Four reference frames are used to characterize the system’s operation. 

The inertial frame, I, is fixed with respect to inertial space. The orbital 
frame, 0, is defined to have axes which are local vertical, orbit tangential, 
and normal to the orbital plane. The computational frame, C, is defined in 
terms of the orbital frame and represents the desired location of the vehicle's 
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body frame, B. Ideally, Che B and C frames should be coincident. Attitude 
misalignment is measured by the small angle misalignment vector, S', between 
B and C frames (4' is the three-vector about which a rotation of jf] radians 
by computational frame C will bring the B and C frames into coincidence). 

Q 

In Fig. 23 a commanded rate to the computational frame is generated, 

(for an earth pointing mission this will be earth rate). This command is sent 
to the stabilization loop and held. The star data, while the rate is held 
by Che stabilization loop, are transmitted to the ground via a radio link, and sent 
to the star-track preprocessor. The star-track preprocessor collects data 
from a known star for one stabilization loop sample time (on the order of 30 
seconds) and estimates the position s"(t) of that star in the star tracker at time 
t. The star-position generator also computes the nominal position for the same 
star, assuming perfect alignment of the B and C frames. 

Any difference between the two star positions is a measure of the 
misalignment between the two frames. This information is fed to a Kalman filter, 
which estimates the misalignment between the two frames, 'l'(t), and at the same 
time, updates the values of any observable system biases. 

The misalignment information serves a dual purpose. It is sent to the 
ground ralselon processor to update current LOS information, and it is also used, 
when desired, to provide a correction rate command, which is added to the nominal 
command transmitted to the spacecraft. The corrective rate will drive 
the current misalignment to zero in a fixed amount of time. 

There are three categories of errors which, if uncompensated, will cause 
the LOS to diverge; they are all of the bias or slowly varying parameter variety; 
errors in the gyro package, such as bias, scale factors, and input -axis 
misalignment; errors in the sensor system alignment; and errors in orbit deter- 
mination and star position. These can all be Included in the state of the Kalman 
estimator. The more-stable errors can be estimated Infrequently in a batch 
mode; the more unstable ones, such as gyro bias, are estimated on-line along 
with the misalignment vector, 4'. The results of a detailed covariance analysis 
of the star-track preprocessor and Kalman filter estimator for a six-state 
implementation that predicts the achievable level of performance are described 
below. 
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Fig. 23 Attitude Determination and Control 
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The Jitter portion of the LOS has Its source In the onboard stabilization 
portion of the system. Tlie main noise sources are the gyro output noise, 
stabilization loop electronic noise, and reaction wheel torque nolae. These 
noise sources must be suppressed through the wide band stabilization loop control 
system and the dynamics of the spacecraft. 

A covariance analysis was performed vrlth all of the above noise sources 
as the excitation. In this analysis, we took advantage of the fact that the 
shape of the output time history from a star detector depends only on the path 
of the star point-spread function as It traverses the star sensor detectors, 
since the Intensity remains constant throughout the traversal. In the analysis, 
we assumed that the star track Is completely characterized by a straight line 
at a known velocity with known point-spread function. Hence, an initial x^, z^ 
coordinate pair and a reference angle measured from a line defined by the locus 
of the star completely determine the output signal. Assuming the electrical 
characteristics are known, and the point-spread function Is an azlmuthally 
sjrmmetrlc gausslan function, the estimate of the Initial star position (and 
angle) Is a recursive update using successive sensor outputs. The covariance 
analysis results tend to be Insensitive to slight modeling Inaccuracies. 

Given a single voltage measurement at specified location on the sensor (a 
sensor node), (v^)> from the previous a priori knowledge we know the least-squares 
linearized solution for the corrections of the track parameters to be given by 

Az f» GhV^(AV) (55) 

' A0 ‘ 

T 

where H Is a 3 by n matrix given by the transpose of the n by 3 measinrement 
sensitivity matrix, H, 

R Is the n by n covariance matrix of n voltage measurements, assumed 
diagonal and proportional to the Identity matrix, and 

(AV) is a column vector composed of the n measurement differences, 
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where C is also the covariance matrix for the estimate of the three-component 
state vector. The square roots of the diagonal elements of C, then, give the 
uncertainties of the x , z and 6 parameters of the star track. 

It can be shown that when each additional node measurement Is combined with 
the current estimate to obtain a better state estimate, the covariance matrix 
relation Is given by the recursive equation 

C‘^1 - C"^ + (57) 

where is the inverse of the previous covariance matrix obtained before the 
new measurements are used. 

Computer analyses have been performed to determine the dependence 
of the covariance matrix on the relative star-track geometry. As might be 
expected, the worst uncertainty in position occurs In the z^ estlimite, the 
direction nearly perpendicular to the star track. If the star Image continues 
to traverse successive sensor nodes with near-zero Inclination angle, the z^ 
coordinate becomes nearly unobservable. In this geometry, all of the output 
signals will be nearly the same for all paths slightly displaced In z from the 
nominal path. By tilting the star tracker slightly this unobservability quickly 
disappears as caa be seen In Fig. 24. 

For a relative star tracker angle of 7 deg, the final standard deviation In 

the estimate of z for a 30-second star-track traversal Is on the order of 
o 

0.05 mil, and somewhat smaller for x^. This result assumes a peak slgnal-to- 
rms-nclse ratio of 10:1. Hence, It seems that the nominal 30-second-duratlon 
collections of data will provide sufficient high-fidelity pseudo measurements 
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to the Kalman filter processing stage for the more time consuming attitude 
determlnatlon/blas parameter estimation function. 


The function of the processing segment described above Is to determine a 
measurement th?t can be used to precisely determine the attitude of the spacecraft 
and provide the uecessary Information for any subsequent attitude control correc- 
tions of unwantec .rifts. A natural byproduct, which Is necessary to achieve 
hlgh-preclslon results, will also be the estimation ("learning”) of at least 
three different types of system biases. This set of parameters consists of 
(1) the three-gyro drift rate biases, (2) the misalignments between the rate 
gyro package a>ces and the nominal body coordinate system axes (each axis 
misalignment being characterized by two direction cosines) and (3) the three 
torquer biases, which will effectively characterize the biases In the onboard 
rate stabilization segment. 

B 

If we define the three vector, w^g, to be the angular velocity of the body 
frame with respect to the Inertial frame, expressed in the body-frame 
coordinates (hence, the superscript), we can write the continuous state equation 
expressing how the misalignment vector, ¥ behaves with time. Before writing the 
equation. It should be emphasized that the misalignments are inf Intesslmal . 

Hence, the state equation will always have at least one error source, even 
If the system model were pe^'^fect. Furthermore, It Is assumed that the differences 
between the actual and nominal rates remain small. The state equation, then, 
for the misalignment vector becomes 

’ 5 ' * + WjB “ ‘^ 1(3 + noise (58) 


where F Is the skew-symmetric matrix 
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B 

When the relation between and the commanded (nominal) rates to the 
onboard gyro package is modeled with unknown biases, the expanded state 
equations can be shovm to assume the following form: 
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5* ■ P(o)*?_) V + B -f r y + r 8 + noise 
1.V Y ® 

B - 0 + noise 

( 60 ) 

e 

Y ■ 0 + noise 
s ■ 0 -f noise 


wlisre B Is the three- vector of gyro biases » y Is a six-component state vector 

expressing the gyro misalignments » and s Is a three-vector expressing the torquer 

biases, and P , P are the ^fect of these on the mlsaligumeat vector, 
y 8 

To construct the Kalmun filter. It Is also necessary to define the relation- 
ship between the observable star position pseudomeausrements and the state to be 
estimated. This equation can be given by 


AX - MT + n 


where AX is the error in the pseudomeasurement two-vector, and Is due only to the 
first three components of the state vector, l.e., the misalignment vector, T, 
and n is the pseudomeasurement noise whose covariance Is provided by least squares 
Inverse above. The measurement sensitivity matrix, M, is given by 


M - 
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where x^, Is the star pseudomeasurement position and f Is the effective focal 
length of the star sensor. 

We performed a preliminary covariance analysis, assuming that only the gyro 
bias three-vector is to be estimated with the Kalman filter. For the six-component 
state vector, consisting of the misalignment vector and the gyro biases, the 
discrete solution for the covariance propagation is given by 
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Here <p Is the transition matrix given by: 
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( 63 ) 


and At is the filter cycle time, nominally 30 seconds. Note that i(i does not change 

c 

with each update since is assumed to be held at some constant value during 
all mission maneuvers. Returning to the above equations, Q is the six by six 
covariance matrix expressing the plant uncertainties, is the previous updated 
covariance matrix, is the extrapolated covariance matrix before updating 

with the current pseudomeasurements, and P is the current updated covariance 

ll+X 

matrix after accounting for the current pseudomeasurements, whose covariances are 
given by R. The matrix is the same function of the current pseudomeasurements 
defined above. 

1 

The main features of the results of our preliminary analysis are shown in 
Fig. 25. For the analysis, a gyro drift rate of 0.01-deg-per-hour-per-day 
was assumed, not an unreasonably demanding requirement. Using only a single star 
in the FOV of the star tracker at any given time, the upper graph of Fig. 25 
shows an excellent transient response for the two sensitive coordinate axes 
misalignments. Tlie third axis will require about 6-hours before comaprable 
performance is achieved. In the lower graph of Fig. 25, the drift rate error 
is seen to achieve 0.01-deg-per-hour within an initial 15-minute time interval. 
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III. THEORETICAL QUESTIONS FOR CONTROL OF LARGE SPACE STRUCTURES 




1. ORDER REDUCTION PROBLEMS 

The fact that the differential equation that result from finite element 
modeling represents undamped harmonic oscillators means that Eq. (7) Is never 
exactly satisfied at any time.* This Is quite distinct from the normal singular 
perturbation approach which we formally used above where the reduced solution 
matches the actual solution to the full set of differential equations for times t 
which are away from the boundary. The reason this convergence is not valid for 
the finite element method Is that there Is no damping. A simple 2 mass model can 
be used to Illustrate this point. Figure 26 shows the model for a simplified 
2 mass example. With the parameters given, the "finite element" model becomes 



The result of transforming via the Cholesky factor of M (since M is diagonal 



“ The matrix $ In Eq. (7) is a transformation that supposedly gives all of the 

components of x In terms of the modes q. When q is truncated, the harmonic 
I motion of the ignored modes never damps, so x is never correct. 
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and if 



( 66 ) 


If the second state (the fast state) is eliminated from Eq. ( 66 ) then the reduced 
model Is 


til “ *^1 ^1 ^2 

« .0032 f^ - .00128 

Assuming fj^ = f 2 = 0 and q^^ 2 ( 0 ) the only Initial conditions on the mode state 
(6j(0) = 62 ^®^ “ then since x « L ^<J >3 the solution to Eq. (67) gives Xj^(t) as 

Xj^(t) = ,06 q^( 0 ) cos t ( 68 ) 


whereas the solution to Eq. ( 66 ) gives Xj^(t) as 

Xj^(t) =» .06 q|^(0) cos /5t + .08 q 2 (^) aos 5t (69) 

If x^(0) = 1 then Eq. ( 68 ) and Eq. (69) can be compared. Figure 27 shows such 
a comparison, and as can be seen the reduced solution Eq. ( 6 S) never converges 
to the actual solution Eq. (69). 

The solution to this problem Is to formulate the order reduction as a "weak" 
convergence. The best way to achieve this is to use the control performance 
measure as the criteria. Thus, if the performance measure for the original 
system is 



then an approximate J is obtained by using a transformation based on diagonal 1- 
zatlon of the cost matrix associated with Eq. (70). The benefit of such a weak 
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order reduction Is two fold: 


1. Periodicity In discarded inodes may be handled. 

2. The dlosed loop dynamics are used as a criteria for 
reduction In order. 

The first use of this form of order reduction was by Meier In Ref. 10. 

2. ORDER REDUCTION IN THE WEAK SENSE 

The procedure we have developed for reducing the order of a dynamic system 
accounts for all of the objections that were raised in Section III-l. Firstly, 
the reduction in order utilizes the performance Index as a criteria tor determin- 
ing the order reduction so that periodic motions In any particular discarded state 
do not cause problems In the singular perturbation. Second, the procedure uses the 
closed loop dynamics. Thus if a high gain system is required the order reduction 
will account for that fact by retaining the "faster” modes In the open loop 
dynamics that are normally discarded. This second step is a significant 
departure from the order reduction described In Refs. 11 and 12. 

The procedure for order reduction differs from singular perturbation 
techniques based on the dynamics of the system where the highest derivatives 
of the fastest states are permitted to go to zero (Refs. 2, 13 and 14). In 
that case, convergence of the reduced model solution to the solution of the full 
order system occurs in an absolute sense. In our case the quadratic performance 
measure Is used to determine the order reduction so that convergence Is In the 
"weak sense". 

To develop the reduction procedure, we start with the open loop dynamics 
(l.e., u ” 0) 

X = A X ; x(0) = Xq (71) 

and ^ 

J = / x"^ Q X dt (72) 

t “ 
o 


since x(t) ■ x(t^) 

J - /[ Q $(t-t^) dt x(tjj) (73) 

o 

or 

J - P(T-t ) X 
-o o -o 

where P(T-t^) satisfies the differential eqiiation 

+ PA + Q (7't) 

at 

which is simply obtained from differentiating Eq. (73) and using the fact that 
<I> and A commute. If T ®, then the steady state solution, P^, of Eq. (74) 
satisfies 

A^P^ + Po„A + Q - 0 (75) 

The solution to the equation (75) is best developed by diagonalizing A 
(assuming A can be diagonalized). Thus, if T is the matrix of right eigenvectors, 
then 

T“^ A T « A 

T 

and if we premultiply Eq, (75) by T and postmultlply by T, then 

tTat“TtTp X + T^’P T T“^AT + T^QT = 0 
00 00 

or 

AT^P T + T^P TA + T^QT = 0 

00 CO ^ 

^ rp T 

and if P “ T I^T and Q = T QT, then P is easily seen to be given by 
_X ~ -1 

Thus, P^ = T P T can be determined by using the eigenvalues and eigen- 
vectors of the matrix A. 
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The matrix Is symmetric, hence It Is diagonalized by an orthonormal matrix 
S , 1 • e • , 

. 0 * 

• ® T 

- SP (76) 

• ® 

n 

where 6, > > . . . > 6 

1 ” ^ “ ~ n 

The order reduction that la now proposed Is based on the eigenvalues of the 
matrix developed In Eq. (76). The reasoning Is that the eigenvalues of Eq. (76) 
specify the amount of control that Is Important as each of theioitlal conditions 
x^ are perturbed, and are therefore a good mi^sure of the effect of the control 
In terms of the desired performance. Obviously the above derivation depends only 
on the initial conditions, since the control was assumed to be zero. When the 
control Is non zero, the same analysis may be used If one assumes that the 
control Is linear. 

X " (A + B K) X ; x(0) - x^ (77) 


A - 


1 

0 6 


0 . 
2 • 


0 0 . 


The last piece we need before the order reduction algorithm is derived 
Is the fact that the closed loop dynamics, when the gain Is known, can be put 
into the form of Eq. (71) and Eq. (72). 

As the first step in the order reduction procedure, assume that the desired 

performance has been translated into the quadratic perforamnce measure 

00 

J = / (x^Qx + u’^Ru) dt (78) 

o 

and that the control is knowri and given by u « Kx. Then Eq. (78) becomes 

00 

J = / x^ CQ + R K) X dt (79) 

o 
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If K la a stable gain In the full state system (x Is a dimension n, which Is large, 
and we desire to find a reduced order model), then the goal Is to find a smaller 
dimension approximation such that the value of Eq. (79) Is approximately the 
same for both the full state and reduced state system. 

Following Eq. (75), the value of the cost matrix P for Eq. (79) Is given by 

_T~ -1 

T PT where 


Pjj - + Aj)J (80) 

and where 

Q - T^(Q + K^RK)T 

T is the matrix which diagonalizes A -f BK 

are the closed loop eigenvalues 1“ 1, ..., n. 

This matrix Is symmetric and hence can be diagonalized as In Eq. (76) so 

P^ - S A (81) 

where 

T T 

SS - S S - I 

In Eq. (81) let us partition P^ in blocks of dimension m and n-m as follows: 


P^ e 

^11 

P 1 
*^12 


CM 

" 22 . 


where P 
P 
P 


11 

12 

22 


m 


Is m X m 
Is m X n-m 
Is n-m X n-m 

is the order of the reduced model which is specified based on the 
error introduced by the order reduction. 
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Thue« from Eq. (81) we get 


^11 ®12 -11 **12 ^11 ® ®11 ®12 
-?21 ® 22 r -^12 ** 22 - Lo ^ 22 J *-^12 ® 22 “ 


T T T T 

or S,,P,, + S" P‘o " A,,S, 


s:,P,, + S‘ P,^ - 


11 11 21 12 “iril ’ “iri2 ^ “21"22 ‘ li '2i 


end 8^2^12 ■*■ ®22^22 " ^22®22* 


T - + cT pT . A 

*^ 12^11 ® 22*^12 ^ 22^12 


By the way In which the elgnevalues 6^ were ordered In A, the elgnevalues in ^22 
are small compared to those in If these eigenvalues are of order e then from 

Eq. (82) as E 0 (the mafaltude of e that is "considered small" in fact will 


determine m) 


gl p ^ p 
'*12^2 ^ ^22^22 


* -T T 
^22 "“®22®12^12 


(83a) 




' 12 - -'’ uSl 2®22 


(83b) 


which gives 


^22 " ®22^12^11®12®22 


since P^^j is symmetric. 

The cost matrix P, from Eq. (83) is therefore 


”^11®12®22 


-T T -T T -1 

r®22^l2^11 ®22®12^11®12®22 J 


(83c) 
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-T T 

«22 ®I2 


~T T 
“® 22®12 


When the Initial conditions are Gaussian zero mean with covariance Q^, then 
E{J) - tr {PQ^} 

- tr{P,, Q + Q P->i Q P'l'' Q > 

11 ''Ojj 12 21 ''Ojj 2. ''oj2 

-T T 

Thus, we assume that '^2 ^ “ ®22^12 -x’ ® partition of. the 

state X Into an m dimensional reduced order subset (x^) and an n-m dimensional 
"residual" subset which are now a linear combination of x^* Using this 
definition, 


-T T 

-®2?12il 


-T T 
-*22®12 


so when x Is used In Eq. (79), 


f x^(Q +K^RK)x dt 
o 


r T 

22®12 


(Q +K RK) 


-T T 
"®22®12 


5l dt 


The state variable mo^^l for x^ is given by substituting Eq. (85) into 
Eq. (79), thus if we call the closed loop state matrix F and partition It as: 


(A + BK) 


where Is m x m 


^22 ^ 
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then 


-1 “ ^11 -1 " ^ 12 ® 22 ^ 125 i * 


X. (0) as before 


(87a) 


and 


-2 


^^21 " ^ 22 ® 22 ® 12^-1 


X2(0) 


- ®22®l2 


(87b) 


are the reduced order state variable models. Since we have assumed X2(0) « 
,-T„T 


~®22®12*1^®^ ’ the performance measure becomes 


J - / x^ (Q +K^RK)x dt 


T 

X P X 
-o -o 


“ x^^(O) 



"l 

p 

“• X 

L-S2?Lj 


L ®22®12 J 


5j(0) 


( 88 ) 


which, if e ->■ 0 becomes: 




Comparing Eq. (89) with Eq, (8A), the order reduction Implied by Eq. (85) gives 
the same performance measure as results when the small terms In the cost matrix 
are aero (e -► 0) . 


To develop the series expansion in e for the cost measure requires an 
tvvi,panslon of Eq. (86) as a function of e. The result of such an expansion will 
give an Indication of the error Introduced by the truncation. This was not done 
during this contract and is one of the Issues that should be pursued. 

The reduced model Eq. (87a) where u is not Kx is given by 


-1 " ^^11 " ^ 12 ® 12 ® 22 ^ 5 l ®1 - 


•I = ^ I _T T 

° '“® 22®12 


[Q] 


T„r I 2l + y ^ H > dt 


-® 22®12 


(90) 
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which gives the reduced order optlroal control problem. The control u Is 
given by u - where is the reduced order gain given by 


and 


0 • 


- r‘^b]’p 

• 

f 

u = K 

1 « 

i - 

- r 

^11 " ^12 

-T T 
**22®12^ 

p. 

+ 

8 

P. 

- 



1 

T 

I 


Q 


..-T..T 


o-T^T 

1^22^12^ 


■®22®12 

mm «i 


P (A 


11 




( 91 ) 


+ P B R"^B^P 


(92) 


The gain Eq. (91) is an m x p matrix, xchlch may be used in Eq. (79) to 
evaluate the performance of the complete n state system as follows: 


E be the n x m matrix given by 

m o ^ 


r I 

m 


0 

L n-m X mj 


T 

x^ = E^ X so Eq. (79), with the gain K^. becomes 


J = / {Q + E k'^ R K e'^} X dt (93) 

- ‘ m 1 r m - 

o 

Obviously, all of the steps leading to the reduced model Eq. (87) can be 

T 

performed using Eq. (93) instead of Eq. (79) with u=K^Xj^=K^E^ x. There is no 

reason, on the first iteration, that the resulting state matrix F in Eq. (87) will 

be the same. Thus, we must Iterate on the gain K and the reduced model 
-1 ^ 

(F, , -F, tS, „S»„) until the result of two Iterations gives the same F. When this 
occurs, the reduced model will have been determined along with the optimal gain 
that accompanies the dynantlc model. 


Let 


then 


76 



Figure 28 shows the iteration that Is described above for the case where 
M X -f K X ■ f are the d3mamlcs (corresponding to the finite element model). 
The matrix A Is then (formally) given by 




0 

I " 


X 


f" 1 

k 

m 

-m“^k 

0 


• 

X 

+ 



(94) 


3. MODE SPILLOVER 

Mode spillover of the control (Ref. 15) refers to the problem that when an 
actuator Is placed on a physical structure. It excites all of the structural modes. 
This can be minimized If the actuator Is constrained to move slowly relative to 
some defined high frequency. The problem is how to limit the control bandwidth 
without adversely affecting the response. 


Because an optimal control design has a l/o) characteristic (see Ref. 7) for 
large frequencies, the fundamental theorem due to Bode Is satisfied (the 
requirement of at least 20 db/decade roll'-off at the zero decibel crossover). 
However this 20 db roll-off continues to the higher modal frequencies where an 
additional l/o) roll-off in each control channel for high frequencies would 
be desirable. To achieve an additional 20 db roll-off, an observer fs Introduced 
Into each control channel. We assume the state to be observed Is the feedback 
gains times the state SljJ* Thus, z the observer state will become the 

control 


z K 


ar 

3iJ 


z of dim. p 


(95) 


The observer equation Is (where the state variable model here comes from Eq. (42)) 


ai* 


r 0 


0 

Li 


(96) 


z = Fz + G I t + K 

- - p- 

where from the observer constraint equation (the matrix K is the observed linear 
combination of states) 


G ■= K 

and L is defined in Eq. (17). 


0 


2 T T -1 

0 0 0 


^FK 
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FINITE ELEMENT MODEL: M£+Kx=Bf m do.l. 
INITIALIZE ITERATION: n = o 


fo. • 

[-M*' i 


; = a stable initial gain 1 = Q ° ^ ^ 


© 

FORM: 


SOLUTION: 


DIAGONALIZE P^; 


A „ P„ + P„ A„ + Q„ = 0 
n n n n n 


" ' \y^\ Vr’A-i-A 


jT . A, O . . O 

o *2 • • • O : . . .> 6 , 


0 PARTITION 


pit S,2l 

L^21 ®22J 


O O . , ,1., 


INTOS,, WHICH ISk X k 
k << 2m 


0 REDUCED ORDER |x,1 ; $2, 2, =■ Sjj *2 ^2 “ *22 ^ ®21 ^1 

MODEL; |s 2 j 

A -= r ° *k ] A r ° 

[m„’<K„-K,2 8,2822’! 

0 DERIVE OPTIMAL GAIN K FOR ii '= ^2l +^S'; 


LETQ„,,=Q.E^K/RK„Em AND A,„ <A - BK„ E J) 


" ['c ] 


^ ® n'==n+ 1, ITERATE ON REDUCED MODELS GAINS 

03S7-028W 


Fig. 28 Order Reduction Algoi ithm 
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The selection of the ’’poles'*^ of the observer (eigenvalues of F) Is made 
In such a way that the Bode criteria Is satisfied at the frequency which Is the 
highest frequency retained in the reduced order model and for frequencies beyond 
that, an additional pole Is introduced. The resulting control system appears 
as shown In Fig. 29. The same transfommtlon used In Section II. 6 Is used here 

In terms of m coordinates of the structure 

4. COMMAND GENERATORS TO STEWING AND SHAPE CONTROL 

The problem of muving a system from one orientation to another may be 
formulated as a control problem as follows: 

"Given an Initial condition x and a terminal condition x«. determine 

“O -r 

the control u such that some performance measure Is minimized (such 

as minimum time, etc.) that also causes the system x « + Du to 

move from x to x,." 

-o -r 

Most of the problems formulated this way lead to nonlinear control laws which are 
difficult to compute - particularly when the dimension of x is large as It Is In 
the structural control problem. Thus, we have formulated the slew problem In a 
different way - an approach that allows the use of linear optimal control 
techniques. 

The assumption Is made that the command Is the output of a linear system 
with an Input that is a step function. Since It Is Important that slew commands 
be smooth, the command generator should be such that as many derivatives as 
possible are zero at t 0 and that near the terminal time of the command 
(t “ t^) the derivatives of the command should all be assymptolcally approaching 
zero. Thus, the command time duration t^ may be selected based on the degree 
of steady state performance desired (how much motion Is tolerable at i:^) . A 
linear system which achieves this level of performance Is the approximation to 
an ideal delay netvrork (Ref. 16) \djose transfer function Is e~ where T Is 
the delay time. The transfer function of the approximate system of order N 
(even) Is given by 

K.) = — ^2 

^N® + • • • + a^'s + a^ 


to provide the modal measurements! ^1 


and their rates. 
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Fiig. 29 Control System with Observer to Reduce Spillover of Control 
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where 


a 


1 



and 


J I Jt 


k/ (j-k)!k! 


By using a partial fraction expansion on T(s)» the transfer function becomes 


T(s) 



+ 



+ 


+ ^N/2 

® -Pn/2 


. ^N/2 

* 

® ^N/2 


(98) 


where denotes complex conjugate, and the state variable form of Eq. (?8 ) Is 
trivially given by (x^ Is an N vector) 


“pj^ 0 0 ... 0 


■*-1 1 

A 


* 

0 p- 0 ... 0 

X + 

b-. 


~c 

1 

• 

• * • 

„ * 


• 

• * 

0 0 ... Pjg/2_ 


- *^N/2- 




(99) 


The roots of the denominator polynomial In T(s) (pj^, ... , Pjj/2 their 
conjugates) are not simply described. They are all complex when N Is even and 
they are not clustered along any path In the complex plane (unlike the Butterworth 
and Tchebyshev filter poles). We rely on the control design calculation to 

formulate Eq. (99) and calculate p^, 1=1 N/2. The step response of 

Eq. (99) (u = 1) Is computed by the addition of an N + 1st state variable In 
Eq. (99). This variable becomes u and since a step has a derivative that Is 
zero, the new form of Eq. (99) Is 


Pi 

0 

0 

* 

Pi 

0 

0 

... 0 0 bj^ 

... 0 0 bj 

Xj, ; y = [11 • • • 10 

0 

0 

0 

... Pjj/2 0 ^2 

(100) 

0 

0 

0 

* * 

• • • 0 Pn/2\/2 


0 

0 

0 

... 0 0 0 

mm 
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T 

with Xj,(0) ■ [0, 0, . . „ , 0 , 1 J . 

Our consputer codes do not use complex arithmetic, so a transformation P 
Is applied to Eq. (100) to make all of the elements real, where P Is given by 


1 1 ' 

-10 j -1 0 

• • » 1 

-10 1-10 

_ _ J.. 

0 

0 


o 



• 

• 

• 

— 
•H -H 
' 

•-I tH 

. 

o o 

■U-^ J^l_ uj 

• 99 

...J 

1 1 

i 1 

• •• i «•» I It* 

j 1 

■ . 

... 


! i 1 J 

n ! 1 

0 j ... j 

111 
1 1 

o o 

-1 

. 0 0 

1 

0 0 j . . . j 0 0 

. J 


It Is readily verified that P Is given by 


n 1 

p o ^ 


1 1 1 
1 
1 

1 -1 1 

1 - 

1 

0 

1 

0 

i j 

! 0 0 
J j 

1 

1 

1 

1 

1 ! 

t 1 

0 i 
1 
1 

1 

-1 

! 0 
1 1 




■"-“t 

-• 

• 

• • 



In the new coordinate system (through the transformation P) , the state equation 
becomes (where In this equation p. = a, + ig. ; j = 1, ... , N/2) 

J J A J 


X 

-c 


“i ®i| 

(a2“0ij^) 

^2 " * 
1 
1 

1 

-h “i j 


1 

0 j ... 
- 1 . 

1 1 

9 

• 

» 

L J 

1 1 

1 1 

^ • 

1 

* * * 

j 

o 

1 

1 

1 

1 

1 

1 

O 1 

1 

! 

■"" ! 

1 

1 ... 

1 

1 


1 

1 

0 0 T 

0 

0 i ... 


!^“ n /2 "* N /2 


) 3 


N/2 

0 


--I-. 

» 

a. 


I 


*N/2 


-3 


N/2 


"N/2 

“n/2 

0"‘ 


K/2 ^ 


I 

j"! 


•k 

21 


^N/2~^N/2 

21 

^N/2~*‘^N/2 
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and the output Is 


y«[l 0 0... 0] 


( 101 ) 


A typical set of plots of Xj^, X 2 » • • • when N ■ 10 Is shown In Fig. 4. 


The next step "‘sing Eq. (101) to design the control system is the deter- 
mination of what In the structure is to be commanded. If the rigid body coordinates 
are the only coordinates to be rotated, then the structural model In the rigid 
body modes are commanded to follow y In Eq. (101). This Is achieved by 
applying a torque that Is the double Integral of y since the rigid body motion 
In modal coordinates Is I 9 T. If the rigid body rotational modes In the vector 
^ (as introduced In Eq. (20)) are the last three components of then the actual 
motion of the physical degrees of ft^eedom are given by 



L $ 



y(t) 


three rigid body coordinates 
are only non zero terms. 


where! 

Xjj(t) are the desired motions of physical coordinates in the finite 
element model (x(t) is the original finite element model 
coordinates) 

L and <I> are as In Eq. (20) 

y(t) is the output of the command generators above. 


The optimal control problem is now modified to cause the actual nodes x 
to follow Xp as follows 


J = / {(x-Xp)^ dt 

o 


( 102 ) 


where 


= 

I 

0 
1 
1 

LIJ 


U 0 ... 0] 


X 

-c 
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and 




^K + Kx - B u ■ j 

As can be seen» the control solution will consist of a gain matrix multiplying 
the states and x. Thus the actuator command will be a linear combination | 

of feedbacks from x and feed forwards from the coimnand generator x^> 

i 

5. ON ORBIT TESTING FOR LARGE SPACE STRUCTURES ‘ 

One of the difficulties encountered when controlling large space structures 
is that the structural dynamics are not known well enough. If a finite element model 
of the structure is used to develop the control system, then the typical errors 
in the mode frequencies and mode shapes that result will lead to a control system ‘ 

that at best does not perform as well as possible and at worst could be unstable. 

The dilemma is that large space structures cannot be built and tested on the 
ground — one must wait until the structure is designed and built in orbit before 
reasonable testing may begin. The problems of this "on-orbit" dynamic testing is 
discussed here in terms of a phase locked loop adaptive spectrum analyzer that 
could provide mode frequencies and mode shapes for control design during the 
initial orbital operations of a large space structure. 

The determination of the structural parameters of a large space structure 
is almost impossible using ground testing. The influences of gravity, aerodynamic 
forces and the difficulty in establishing the thermal gradient that will be 
encountered in orbit all contribute to uncertainties (even when these effects 
are analytically extracted from the test data). Added to these problems 
is the Impossibility of even assembling the actual structure on the ground so 
that it maintains its Intended shape. Thus it becomes Important to consider 
testing the structure when it is in orbit. 

The use of modern control theory to develop control systems for large 
structures requires that a detailed model of the structure exist. Due to the | 

difficulty with testing the structure on the ground, the control design must be 
deferred until the system is in orbit and the structural model becomes available. 

Use of a digital control system is then ideal, because a rigid body low band- 
width controller may be used initially (during d 3 mamic testing) . After the control 
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algorithm Is developed using the test data* the original low bandwidth controller 
may be easily stripped out of memory and replaced with the control system designed 
using the structural data. Thus» the question becomes one of developing a 
structural model. 

Since the structural frequencies are discrete, a Fotirier transform of the 
output of a sensor mounted on the structure will show discrete lines at the modal 
frequencies. The use of a spectrum analyzer is thus a possible way of developing 
a structural model. The problem with this approach is that a decision process 
must be appended to the output of the spectrum analyzer to allow the discrete 
mode frequencies to be selected. This decision process must be automated since 
one would expect to do the dynamic testing periodically to update structural 
mode data as the structure's properties change. 

A natural device for automatically selecting the mode frequencies is a 
Phase Locked Loop (PLL) . This device suffers from some problems when the 
frequencies of two modes are close together because it will alternately pick 
out one or the other frequency. In addition, the loops are sensitive to certain 
noise processes. The characteristics of a PLL are described and a method is 
Introduced that will allow the loop to be better tuned to the characteristics of 
structural dynamics. The theory for identification of structural parameters that 
are built Into the new phase locked loop, and some results from a simple computer 
simulation of the loop are also shown below. 

The technique we are describing here uses an optimal filter in the loop 
to tie the loop operation to the known characterltcs of the structural dynamics. 
Reference 1 is the closest application of such an approach in the literature. 
There, a Weiner filter is developed to give an optimum filter for a single 
sinusoid, whereas here a Kalman filter with a frequency identifier is used. 

Figure 30 shows a phase lock loop as it is uorraally configured. The 
operation of the loop relies on the fact that the result of multiplying two 
sinusoids is sinusoids at the sum and difference frequency of the two sine waves. 

Thus if the input sinusoid is sin (oi^t + <^) and the output of the voltage control- 
led oscillator (VCO) is sin -f 0) the input to the lov/ pass filter is given by 

E(t) = sin t+<(») cos t+0) = 

o o 

h sin((a)^-O)^)t+<|>-0) + h sin( (w^+oi^) t+(>)+0) (103) 
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Fig. 30 Phase Lock Loop Block Diagram 



Fig. 31 Modification to Phase Lock Loop to Allow 
Estimation of Mode Frequencies 
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and the output of the low pass filter will be the difference frequency component 
of Eq. (103). The VCO will change frequency following the low frequency term 
until sin { (u)j^-u^)t+((i-9 }■ 0 at which point the loop is locked where and 

• B. The assumption that the input signal is a single sinusoid is crucial to 
the op-^ratlon of the phase locked loop. If a second sinusoid with a frequency close 
to the primary sinusoid exists in the input then the VCO will not stay locked 
on the primary frequency and will alternate between the two frequencies in a random 
way. 

Figure 31 shows the modification to the basic phase lock loop that was 
first shown in Ref. 9. The fundamental feature of this loop is that the optimal 
filter will have a variable bandwidth vrlth multiple notches. The variable 

A 

bandwidth is a consequence of the convergence of the loop frequency, to the 
frequency contained in the signal u)^. The multiple notches in the estimator comes 
from the apriori assumption of multiple frequencies in the signal. These 
t\‘io properties overcome the major objecfcionft to the use of a PLL for structural 
frequency determination. 

The problem of identification of the coefficients of a linear differential 
equation may be formulated as a nonlinear filtering problem since the unknovm con- 
stants may be assumed to satisfy a differential equation where the constant's 
derivative are zero (Refs. 18 and 19). 

For the phase lock loop used to identify the structural mode parameters, 
the assumption is made that the underlying structural systems is modeled by 
a finite element model of the form 

MS + Kx = f (104) 

A series of transformations are used on Eq. (104) to obtain the modal form of Eq. 
(104) as follows 

• Transform from x to z using the C'holesky factor of the mass 

T 

matrix M. Thus if L is lower triangular and LL = M, then 
T 

defining z = L x and substituting in Eq. (104) gives 

2 * - z + l"^ f (105) 

• Transform Eq. (105) to diagonal form using the orthogonal trans- 

T 

formation g = z to give 
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3 - - g (106) 

2 

where Q Is a diagonal Tnatrlx which has the square of the wode 
frequencies along Its diagonal. 

In Eq. (106), each mode satisfies an Independent differential equation of the form 


d. ■ - 0 ) . q. + b, -f, + ... + b. f (107) 

^1 1 H 11 1 Ip p 

T -1 

where the b . . are the terms In the matrix <t> L . 
ij 

Let us consider a single mode with p = 1 (l.e., only a single force Is 
applied to the system). Since we are postulating a dynamic test mode, we can 
use the actuators on the spacecraft to excite the structure one at a time. The 
unknowns In Eq. (107) then are and . Since the solution to Eq. (107) 

Is of the form of the sum of sinusoids multiplied by the modal Initial condition 
plus the forcing term, the b^^ Is a scale factor on the amplitude of the steady 
state oscillation indui-ed by the force f , . By adjusting the amplitude of f_, 
the coefficient b^j^ may be made unity. Thus, the only unknown parameter is 

If we write the equation (107) in state variable form, where we further 
2 

assume is also a solution of a differential equation we get 




where in Eq. (108) the white noise terms w^^ and W 2 represent the uncertainty 
In the problem. Wj^ is the vibration noise that is exciting the structure and W 2 
is used to change the rate of convergence of the estimator (o^^ and 02 are the 
standard deviations of these noises and is assumed known). 
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2 

In Eq. (108), since o)^ Is a state (the 3rd component of the state vector), 

the equation Is nonlinear. To linearize Eq. (108) and estimate the coefficient 
2 

CD^, an approximate nonlinear filter described In Ref. 19 Is used. This filter, 
closely related to the extended Kalman filter used In inertial navigators uses 
terms up to second order in the Taylor series linearization of Eq. (108). The 
filter uses as a measurement the standard phase lock loop measurement (as In Eq. 
(103)) thus 


where 


yj(t) 


y^(t) “ cos w^t yj(t> 

Is the estimate of the mode frequency 

Is the measurement at the physical degree of freedom 
on the actual structure labeled Xj and Is given by 

Op 

y. (t) = s «<).. q.(t). 

J 1»1 1 


(109) 


Since fj^ was adjusted to give unity mode motion, will be the actual 
amplitude at the frequency contained in the measurement y^(t). This means 
that once the frequency is at the modal frequency u)^, the output of the phase 
loop will be the mode Influence matrix. Thus if a full set of m modes are 

desired, altogether m measurements with m loops at each measurement point are 
required (a total of m phase lock loops connected through a single filter) as 
shown In Fig. 32. 

The resulting nonlinear fllcer is given by 


d 

dt 








p 

*^23 

0 


+ K{cosw^t yj(t) “ cosw^t(()^^qj} 


(liO) 


’Si 
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Fig. 33 Estimate of Mode Frequency for a Single Mode. Effect of Noise Vari- 
ance Assumed for the Measurements (o2) on the Rate of Convergence. 



90 






where P 22 element of the estimation error covariance matrix P and 

K is the "optimal gain" given by 
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- P [ 

COSU)^ 
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0 0] 
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0 

'2 
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- 0 

0 
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2 

where cr^ Is the variance of the measurement noise on the sensor that is measuring 
the physical motion yj(t), 

A single mode system (m » 1 in Eq. (106)) was simulated to determine the 

operation of the optimal loop filter. The result of estimating the loop frequency 

Is shown in Fig. 33. One of the unique characteristics of this system Is the time 

varying dynamics of the loop filter. As a function of time, the loop filter tends 

to start out with a high bandwidth which gradually gets smaller as the loop fre- 

2 

qutncy estimate gets better. In the absence of a noise on the coefficient 
(l.e., when o^ ■> 0) the loop filter bandwidth goes to zero in steady state. 

Thus the parameter 02 can be used to adjust the steady state filter bandwidth. 

In practical application each filter "will have a variable bandwidth and 
will have multiple notches (at the frequency of the modes not being estimated). 
Furthermore, the actual loop frequency, once estimated, will not be changed 
unless there is some reason to believe the estimates are Incorrect. We are 
currently building a large simulation code that will estimate multiple modes to 
test the method on a relatively large problem. 
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6. STRUCTURAL DAMPING IN FINITE ELEMENT MODELING 


The finite element model Is an approximation to the underlying partial 
differential equation of the structure it models. As such the full order 
finite element model is the only vehicle available for verifying the stability and 
performance of the control system design. Since the low order design model 
usually has some dam:’ng assumed, there actually Is a measure of stability Imposed 
on the full order design. This stability Is a consequence of the way the damping 
Is Introduced, and as such It Is Important that the control designer understand this 
modeling. 


In order to preserve the mode transformation developed where the damping Is 
assumed zero, the damping matrix that is used to model the damping on the structural 
degrees of freedom Is assumed to be (Ref. 20) 


where 


C « M 


“ -11 
E a (M K)^ 

1=1 


C Is the damping matrix (see Eq. (3) In Section Ill-l) 
M Is the mass matrix 

K Is the stiffness matrix 


are coefficients which are determined by the amount of damping 
desired on the various modes 
m Is the number of modes In the reduced order model 


(I’*’) 


Using the notation of Section II- 3 the dlagonallzatlon of the matrix C 
uses the same transformation that diagonalized K and M. Thus If 

T T 

2 =■ $ L X 

where 

T 

M = LL 
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Chen 


T -1 -T 
rL C L *4 


£ (M'^K) ... (M"^K)L"^4 

1-1 ‘ ' rTN<-- ' 


1 terms 


® T -1 -T T -1 -T 

I a ((f^L ^KL $) ... (4> L ^KL ^4>) 

i-1 


I a (n ) 
1-1 ^ 


2 . 4 . 2m 

a«Wo + • • • w 
112 2 m 


2.4, 2m 

a,o)_ +a„u)rt + ... +a (1) 

1 2 2 2 m 


0 . .. 0 


... 0 


2 2 2 
0 ... a,(D^+a_u>^ + . . .+ot 0 )^ 
1 m 2 m m 


T T 

Thus, the transformation <I> L diagonalizes M and K and therefore also C. To 
determine a 2 , ... , the levels of damping in each mode Is determined 
(l.e,, the terms 1 = 1, ... , m) 


Then the equations 


^1 ” 


2 , 4 . . 2m 

a,u). + a„u). + ... + a u). 
11 2 1 ml 


2(i) , 


1-1, 2, . , . , m 


are solved for ... , a^. Tills set of equations may be written as 

2m 


2 4 

U)l (D^ 


2 4 

U) 01 

m m 


2m 

i 

m 


Oi 


m 


201^^1 
2(02^ 2 


2u) 5 
m^ m 


(114) 


(115) 


(116) 
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and hence 


Sc 

. . . I-. 
1 

■ v-i 

• 

t 

t 



m m 

U J 


where V is the vandemonde matrix of Eq* (116), 

Ztotice that there is now an implicit damping for all of the higher frequency 
modes, m + 1 to n. These damping coefficients grow with frequency. Since 
the same transformation applies to the full nth order system, the damping of the 
higher frequency modes is also given by ^ 4 . (115) where 1 ■ m+1, ..., n. 

Obviously as the frequencies o) get larger, » which means that these 

higher frequencies are no longer going to be oscillatory. Thus, if this model 
is used to verify stability of the control design, an unrealistic damping is 
actually being used on the modes v;hloh are expected to be problems. 



IV. SIMPLIFIIJD RESIGN EXAMPLES 


The two mass system of Fig. 26 will be used here to describe the linear 
discrete time optimal control techniques that were used to develop the control for 
the OCDA in Section I. This example has most of the features that were needed 
for the understanding of the larger problem, b»;t it is sufficiently simple 
that all of the matrices etc., may be written explicitly. 


The problem is to control the tv«> mass system using the force fj^ only. 
Thus, the "finite element" model is 


"100 

0 “ 


Pk 1 
^1 


"1780 

-1780 * 


■-i' 

+ 

’h' 

0 

343.8 




-1780 

Ml 

4194.4 


.’‘2. 


0 


(117) 


STEP 1 ; COMMAND INPUT 


Let us assume that the desired response of Eq. (117) is such that the position 
of the dominant generalized coordinates follows the profile shown below 



-t/x, 



where Xjj(t) denotes the desired position of both and X 2 « Clearly, with 
only one force alplied to mass 1 , the positions of masses one and two cannot both 
be Increased from their initial values by 1 unit since, in steady state, if 


Xj^ = 1 and 


0 

IL. ^ 


1780 -1780 


-1780 A194.4 




1780 


1780 

4194.4 


•) . The procedure we 


from Eq. (117) gives X 2 = 4 Y 94 ' ' 4 ^1 “ 1780 (1 - 

described in Section III solves for and X 2 using a weighted least squares 
approach so that and X 2 are both moved approximately the desired one unit. 
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The motion x^Ct) is Imparted to the domi.nant mode, that la, from Eq. (117) 


and since 


- - 5q^ + .06 


, “T , -T . T 
X"L z*L ^ 9 


m mm 

0.1 0 


0.6 0.8 


0.06 0.08 




3 ■ 


0 0.054 

mm m 


0.8 -0.6 
mm «■ 


.0432 -.0324 

M 


which gives 


“ 0.06q^ and X 2 * .0324 q^^ 

If qj^(t) “axjj(t), then Eq. (118) gives 

. -t/Tj3 “t/Tj. 

fj(t) « ^ e *"+5(1-6 *")} a/. 06 


a 


-t/t 


1 — t/ T TJ — I Q 

Qg {5 - (5 + “^) e } “ {5-14e } a/. 06 

"d 


(118) 


(119) 


now the amplitude a of the command x^^ can be determined by solving for ot such 
that x^ and X 2 are as close to 1 (in steady state) as possible. Thus we 
want to solve 

5a 


[0.06 .0324] 


.06 


[1 11 


( 120 ) 


which can only be solved in the least squares sense (this is two equations in one 
unknown) to give (following Section III-4) 


5a 

.06 


= [0.06 0.0324] 


.06 + .0324 


1 

LiJ 


^.06)2 + (. 


0324)' 


.0046 
a = .238 


19.872 


( 121 ) 
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Therefore, the desired motion is 

-t/Tj) 

Xp(t) - .238 (1-e ) 

which may be modeled as the differential equations 


k 


D 


1/td Xjj + 


.238 

"d 


with 


s(o) ■ .714 

Xp(o)" 0 


- 3 + 


e 


s ■ 0 


STEP 2: CONTROL AND COMMAND MODEL IN STATE VARIABLE FORM 

The combined slew and reduced order model in state variable form can be 
written in terms of an augmented state variable p defined as 
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0“ 
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^1 
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X 

so that £ =» 
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L J 


where 


'qi(o)" 


"6 x^(0) + 14.81 X2(0)“ 

qi(0) 

= 

6 Xj^(O) + 14.81 * 2 ( 0 ) 

X^(0) 


0 

0 



s(0) 


.714 

U ^ 


T T 

The initial condition on p is obtained from the fact that 3 = ^ L x. 


(122a) 

(122b) 


(123) 
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STEP 3: PERFORMANCE INDEX 

The optimal control is required to minimize the difference between the actual 
response x. (t)» x»(t) and the desired response x. (t) and x» (l). Thus, let 
us assume that we want to minimize 

00 2 2 
J " / - X (t)] + u, [x«(t) - X, (t)} + rf-^(t)} dt (124) 

o D D 


The reason this performance measure is over an infinite time interval is that 

a constant gain system is desired. The only parameters in J are q and r since 

any multiplier on the first term can be factored out of J (q and r are used to 

adjust the relative match of x. and x~ with x^ and x„ and to adjust the maximum 

^ ^ u u 

force applied). Since x. and x» as well as x^ and x,. are modeled by the order 

I ^ *^D 

reduction as functions of q^, q^^ and x^^, the performance measure Eq. (124) may 
be written in terms of £ as follows 


The vector 


Xo 


'D 


can be written in terms of £ using Eq. (118) as 


f* 

Xi 


'■ .06 

0 

0 

0" 


r „ 
‘Ji 

^^2 


.0324 
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1 
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' S 1 
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But the Integrand in Eq. (124) is 
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L% 


0 


q 

-1 0 
0 -q 


-1 o’ 

0 -q 

1 0 






+rf: 
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Hence, in cenus of the reduceit state vector g, J becomes 


.06 .0324 


.06 .0324 -1 0 

0 0 0 -q 


0 -q 


.06 0 0 0 

.0324 0 0 0 

0 0 .06 0 

0 0 .0324 0 


+ rf 


(.06)^+(.0324) q 0 -(.06)^-(.0324) q 


-(.06) -(.0324) q 0 (.0S)^+(.0324)^q 


^ ^0 


•f rf. 


(125) 


This performance measure, when applied to the reduced order system described 
In Eq. (118) will cause the control system to be structured as shown In the 
block diagram on the following page, where K 2 , K^, are the gains that 
are derived to make f^^ ■ via the optimal control derivation. 
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STEP 4: OPTIMAL CONTROL 

The steady state optimal control comes from the Pdccattl equation 
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m 

(-10p^2-Pr-8p^2^> 

• • 

(- 5 P 22 +P 11 - 8 P 12 P 22 ) (2Pj2”BP22^^ 

symmetric 

• • 

<-5p23“3Pi3-Y“3Pi2P23) ^Pi3“2P23“3p22p23^ 

<-6p33+Y-3P32^) 

(-5p24+Pj3-6r'i2P„) Pi«+P23‘®'>22'’24 

m 

(-3P34+P33“6P23P24) (2p3^+P^^^) 


« (127) 

where y ■ (.06)^ + (,0348)^q 



This matrix has been partitioned Into four 2x2 blocks so that the point 
can be ntadci that the optimal control feedback gains are independent of the feed 
forward gains. This Is obvious In Eq. (127) because the 2x2 block In the upper 
left depends only on Pj2* ^22’ fhls block may be determined 

Independent of the other blocks. Since the feedback gains are given by 
K ■ " « B ^nd since B is zero everywhere except In those positions corresponding 
to upper l??ft most 2x2 block of the gain also only depends on the 

el^vjents p^,, Pj^2* P22' Striving Eq. (127) gives the upper 2x2 block and 

the feedback gains as: 


f / 2 ^ 

>12 “ ~2 ~ 'l+(.06^/5r) (. 

. 06 ^ \ 


06^ + .0348^q] 


/lx p 


12 


22 .06 


; p 


11 


/i+ (.( 


06^/5r)“ [.06^ + .034s2q] 


the control gains are 

- 5/. 06 (l - / 1 + (.06^/5r)^ [.06^ + .03482q]j - 5/.06(l-A) 

/2(~K,) 

h ” " - /lO(A-l) /.06 
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Since the nrade we have retained Is the differential equation 


and since f ^ 2^3 ’ ‘ closed loop dynamics are given by 

~ .06 + (5-.06Kj^) qj^ « 0 

or 


+ / lO(A-l) q^^ + 5A qj^ » 0 

Hence the new undamped natural frequency and damping of this mode becomes 

u) “ / 5A 


and 

5 - € 

Since 

A = / 1 + (.06^/5r) t.06^ + .0348^q] , 

the mode Is damped to .707 when q is large relative to r , l.e., whenever the 
control saturation is not important. On the other hand, as r^ gets large 
(less and less control authority is permitted relative to q, the mode is less 
dampled and in the limit no control at all Is exercised. Note that as A -► «, 
the closed lo<Dp undamped natural frequency (uj) Increases as the damping gets 
closer to .707 which is exactly the problem \d.th the order reduction as it was ' 
performed because when A » 5, the mode retained (q^ ) crosses in frequency with 
the discarded mode (q 2 > whose frequency was 25 rad/sec. 

To make these points clear, a root locus plot is shown in Fig. 34 that 
shows the closed loop pole locations of the design as q and r are varied. 

The feed forward gains which determine the optimum response to the command 
such that x^ and X 2 are as close to 1 in t!ie least squares sense are given by 
K 3 “ “*06/r P 22 and = -.06/r P 2 ^. From Eq. (127) P 22 and P 2 ^ are given by; 
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with the feed forward control gains which give a force equal to K v + K s 
is the following block diagram ^ 
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Hence, the transfer function from s to f^^ Is given by 


fo (s) 
S 



K^{S+(K3/K^+3)} 


The last point that must be made Is that In Implementing the feedback control 
derived above, the measurements will be on and X 2 > hence since the feedback 
force Is given by 


fl -K 


- IKj Kjl 

■^l‘ 


-V 


^1 


(128) 


but since ■ 6x^ -f 14.83 X 2 the Implementation of this control seems to require 
four measurements (x^, Xj^, X 2 » and Xj). Actually, the measurements need only 
be Xj^ and x^ as In Eq. (49). 

Via a series of simple block diagram manipulations, the above closed loop 
system becomes 


I 1 
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which has a closed loop root locus that is always stable so that, for this 
problem, no additional stabilization need be done and the problem may be 
considered complete. 



STEP 5: DIGITAL CONTROL 

The state variable model Eq. (123) is used to derive a discrete model. Since 

t 

x(t) = ‘Kt-t^) x(t^) + / 4>(t-t) b f (t) dT 

t 

o 

if t = (k*"l)At and t^ = IcAt, then from Eq. (35) the discrete system becomes: 

At 

-ld-1 " -k ^ - fl(I^t) 
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where 


<t(At) ■ 


I //5 sln/SAt 0 


-/5" sin At cos/SAt 


At 

/ (!• (x)dx b 


.012 (1-cos v^At) 
.027( sin /5At) 


- 


If the sample time At is 1 sec., then 


-.617 .352 

-1.76 -.617 


.0497 .317 


This model is used with the discretized performance index to design the 
optimal digital control. 


V. COKCLUSIONS AND RECOMMENDATIONS 


This report represents one step on a relatively long road that will provide 
the technology for designing control systems for large space structures. This 
technology is of extreme importance for without iti the schemes, ideas and dreams 
for a large number of space missions will never reach fruition. The work here 
has demonstrated that control systems may be designed that provide structural and 
attitude control when the control specifications are not severe. It has also been 
shown that many of the technology items: structural modeling, modeling of damping, 
determination of the structural dynamics in orbit, slew command and control 
spillover, all have an effect on the con; rol solution. One of the significant 
contributions of this effort is the realization that a control system should 
exploit the ability to "play" one disturbance off against another. This synoptic 
design approach can pay very great dividends in the use of actuator fuel and 
performance. Also a synoptic design is best developed using linear optimal 
control techniques because of the natural incorporation of disturbances, dynamics 
of the structures, rigid body torques and command generator dynamics in the design 
model. If there is one feature that we believe is important in this work it is 
that the linear optimal control design philosophy has been used throughout the 
work. The design incorporates the concepts of stability margin, phase margin and 
in general the robustness properties that are usually considered "classical" 
but using the terminology and theory of linear optimal control. 

Many of the steps taken here are tentative. In particular the "on orbit 
dynamics test" procedure must be developed so that multiple mode frequencies 
may be estimated. The order reduction in the weak sense should be attempted on 
a large problem to demonstrate its effectiveness. Finally, the whole question of 
Kalman filter sensitivity vs. the gain margin achieved by the introduction of 
an observer must be resolved. 

In the course of this effort, three technical papers were presented. 

These papers are Refs. 9, 21, and 22. 
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APPENDIX A 


INTRODUCTION OF GRAVITY-GRADIENT AND 
ROTATION TERMS INTO STRUCTURAL MODELS 


As indicated in Figure A-1, the analyzed satellite is represented by a 

number of lump masses. A typical mass m. , has a moment-of-inertla matrix 

( 1 ) ' 

[I^] in the local body axes x Point 0 is the location of the cm of the 

satellite, assuming that it traveled in the nominal orbit. This orbit is 

specified by the user. The user also specifies the initial position of the Z 

axes and their constant angular velocity These axes rotate at the nominal 

angular velocity of the satellite (e.g., the orbital rate for earth-pointing 

satellites) so that motions of the satellite relative to the Z axes are small. 

{a^} is the undeformed location of m^ in the Z axes, and {v^} and are the 

translational and rotational deformation, respectively, of m. as observed in 

( 1 ) ^ 

the Z axes. In the undeformed system, all of the axes are parallel to the 

Z axes. 

The equations of motion are 


[ra]{»} + [k]{x} = {f} 4- {f } - {g} , 

O 


where [m] is the system mass matrix, 


[m] = 






f 


where [j^] is the 3x3 Identity matrix. In Eq. (1), [k] is the system 
stiffness matrix, {x} is the deformation vector, 


(A-1) 


(A-2) 


The tilde under a symbol is used to designate a nonscalar quantity such as 
a vector, a matrix, or a set of Euler angles. 


A-1 




Fig, A-t Representation of Satellite and Coordinates Used 


A-2 



(A-3) 


{x} 



-M 

-n 





{£} contains the forces and torques on the system other than the gravity-gradient 
and stiffness loads. 




(f) 


T . 




{f } contains the gravity-gradient forces and torques, 
R 


(A-4) 




V 

iSi-' 

i ■■). 

F 

-gn 

Tgn-^ 


(A-5) 


and {g} contains the mass times acceleration terms attributable to rotation. 


{g} 


X I 

r . 


rn 

V. -n J 


(A-6) 


A-3 


The gravity forces in Eq. (A-5) are expressed In the Z axes and are 

8 

(Fg^) - - ra^ ^ ({p^}- 3({e^)’^{p^)){e^}). (A-7) 

where g^ is the acceleration of gravity at point 0, and {p^} Is the location of 
10^ In the Z axes; l.e., 

{Pj^} - (v^) + {a^) (A-8) 

The gravity torques In Eq. (A-5) are expressed In the axes and are 

‘V’"’* (A-9) 

where [F( )] Is the cross-product function defined In Appendix and 

{ej} - ((I 3 ] - [r(4^pHe^}) (A-10) 

The components in Eq. (A-6) that im)dify the force equations are expressed xn the 
HI axes and are 

«v 

{g^} - (2[r(J2) H<r^} + [r(fl)l^{p^}) (A-11) 

while the components in Eq. (A ) that modify the torque equations are expressed 
In the x^^^ axes and are 


{h^} - (A-12) 

where Is the angular velocity of m^ in the axes; l.e., 

+ ([i3l-[r(<|)^)l{n} . (A-13) 

Reduction to Modal Form 

The modal matrix is arranged as follows: 

m « [/ f f] , (A-14) 

f t 

where [(j) ] contains the flexible modes, [(j) ] contains the rigid-body translation 
modes, viz. 


A-4 


'.3 


?3 


I*'l 


(A-15) 


where 0^ is the 3x3 zero matrix, 
viz. 


r 1 

] contains the rigid-body rotation modes, 




•r (a^) 


13 


-r (a ) 

~ n 


l3 


With 


(A-16) 


{x} - [0HU 


(A-17) 


the equations of motion In modal form are 


[MHO + [CHC) + [KHO - (O 


(A-18) 


[M] Is the modal-mass matrix, 


[M] - 


m 


I3 


(A-19) 


A- 5 


where the p^'s are the modal masses corresponding to the r flexible modes that 
are used, 


^1 " ^*^1^ ; 1 - 1, , r, (A-20) 

m Is the total system mass, 

m - J; m^, (A-21) 

and [I] is the total moment of inertia of the undeformed vehicle about its cm 


relative tc the Z axes 


[I] - E (U^J - m^[r(a^)]^) . 


[K] is the modal'-stlf fness matrix 


(A- 22) 


[K] “ [(|)]nk][({)l« 


K, 




K 


where 




1, • . . , r 


,th 


<A~23) 


(A-?A) 


with 0 )^ equal to the 1 flexible frequency. [G] is a modal damping matrix that 
has been added to the formulation at this step: 


[C] 


L ?6J 


(A-25) 


{=} is the equivalent modal force which contains the corlolis and centrifugal 
reverse-acceleration forces -{g} as well as the physical forces {f} + {f } ; vlzj 


{ = } = [<|.]"({f} + {f } - {g}) 

Motion of Rigid -Body Axes 

The modal-displacement vector is partitioned as follows! 


(A-26) 


A-6 


(A -2 7) 


where (C). {p)» and (y) are the flexible, rigid-body translation, and rigid-body 
rotation coordinates. It can be shown that tp} locates the cm of the satellite 
and satisfies the relations 

t 

m{p} - ^ (A-21 

As Indicated in figure A-2 , ip} locates the origin of a set of axes known as 
the rigid-body (or mean) axes of the satellite, denoted as the x axes, (y) 
orients these axes relative to the z axes. For a given deformed shape of the 
satellite {x},{y) may be defined by the relation 

T 

[IJ{y) [mHx) (A-2‘ 


The total motion nuiy be decomposed as follows: 


{x} « ff‘^Hc} + 


where 


I'l' KC) Is the nexibie motion relative to the rigid-body axes 
Is the translation of the rigid-body axes 
HyJ is the rotation of the rigid-body axes 
The components of Eq. (A- 30) are 

Cvi) = + {p} - fr(a^)]{Y) 0 

C'l-'j) = + fY? 0 

ft f 

where ) is the partition of [<)' 1 corresponding to the translation of m . , 

fr^ f ^ 

and ] Is the partition of fi}> ] corresponding to the rotation of m^. The 

flexible contributions 1, Eq. (A- 31) and Eq . (A-32) are called {p^} and {6^^); 

l.e. , 



Fig. A-2 Vector Repr^nlation of Deformations 




(A- 33) 


{ 61 ) (A-34) 

{v^}, {p}, (y)* {u^}* and (0^^) are expressed In the Z axes. Figure A- 2 shows 

^ *** 

the components of the motion In vector form, a^ (which locates the undeformed 
position of m^ as seen by an observer fixed In the x axes) Is the vector a. 
rotated through the angle y* 

The last six equations of the set Eq. (A-18) are the overall equations of 
motion for the system, l.e., the equations for {p} and {y}. The first three of 
these equations are 

m{i5) - Zii ) + {F„} -2m[r(fl)Hp} - m{r(J?)]^{p} (A-35) 

where {F } Is the gravity-gradient contribution to the total load; l.e., 

O 

{F } = - mi ({p} - 3 ({e^}’^{p}) {e^}) (A-36) 

“ o 

The last three euqtlons of Eq. (A-18) are 

[I]{Y> = {T^} + {Tg} - I (2m^[r(ap][r(JJ)]{v^} 

+ [r(a^)][r($^)]2{p^} + [I^][r(n)]{t|»^} 

+ [r(a)^)][I^]{(i)^}) (A-37) 

where is the resultant torque of all of the external loads, other than 

gravity loads, about the cm that the system would experience if it were rigid, 

» 

and {T } io the contribution to the equation attributable to gravity, i.e., 

O 

{TJ = [< 1 .'']’^ {f } 

O S 

Total Gravity-Gradient Torque on Satellite 

{T } is an approximation to the total gravity-gradient torque on the satellite 

O 

however. It does not Include certain effects attributable to vibration. The 
total gravity torque can be obtained more accurately from the following equation. 
This torque is expressed in the rigid-body, or X, coordinate system. 


A-9 


.-■m' 


{Tg) - ([| 3 ] - [r(Y)]) r[r(rp]{Fgj}+E ([|H-[r( 0 ^)]){Tg^} (A-39) 


where 


{r^} - {p^} - {p} 


{0^} « - {y} 


(A-40) 

(A-41) 


Gravity Constants 

For reference, note that 


O 


R 




(A-42) 


where is the earth's gravitational constant, is the earth's radius, and gj, 
is the acceleration of gravity at the earth's surface. Equation (A-42) is useful 


for calculating g^ and K^. 


Particularization to Nominal Circular Orbit 


For a circular orbit 


g = u) R 
o 00 


(A-43) 


The components of the unit vector along the local vertical in the 2 axes, 

are needed for use in several of the equations. In the A axes s-own in Figure 

A-3, {e }is called {e } where 
° ° A 


{e } = j cos(o) t + 5) 

° A \ ° 

j sin(o)^ t + 6) 


(A-44) 


Then, {e } in the Z axes is 
o 


{e^} = [T(a)](Tr(6)He^} 


(A-45) 


A-10 




NOTE: THE EULER ANGLES 0 ORIENT THE 

2° AXES RELATIVE TO THE A AXES 



0357-037W 


Fig. A-3 Nominal Coordinates Relative to Earth-Centered Coordinates 
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A-11 

OMI 





where [it ( 6)1 is the coordinate thransformatlon defined In Appendix B for the 

O 

Euler angles B - ( 6 j^, 62 * ^ 3 ^ orienting the Z axes relative to the A axes, 

o * fi t . 

O 

{fi} Is given In the Z axes as 


{Q} - {nl, 

O 

and [T(o)l Is the transformation from the Z axes to the Z axes due to a 
rotation o about an axis, along {fl}, fixed In space; viz. 



■ l-A(l-ni^) 


Jin2n3 - sr)2“ 

[T(o)l “ 


l-d(l-n2^) 

ilTizIs + srij^ 

where 

ariiria + STI 2 


l-Ji(l-n3^) 

s = sin 0 , c 

= cos 0 , A 

= .l-c 



(A-46) 


(A-47) 


(A- 48) 


(A-49) 
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APPENDIX B 


DEFINITIONS OF MATRIX FUNCTIONS 

Cross-product matrix; 


For any vector x * 

[Xi X 2 x^]"^. 


0 -X3 X2 

fr(x)] = 

X3 0 -Xj 


-X 2 X3 0 


Euler-angle coordinate transformation; 

For any set of ordered rotations y = (y^ y^, Y 3 ) about axes 1, 




1 

0 

0 


[a(y)J = 


0 

cosy 

siny^^ 




0 

-slnYj^ 

cosy^ 



r 






COSY2 

0 

-siny 

[B(y>] = 


0 

1 

0 


[_slnY2 

0 

cosy 


” COSY3 

slny^ 

0 " 

[C(y)] = 

- 

■siny 

j cosy 3 

0 



0 

0 

3 

J 


The total transformation is 

[tt(y)] = [C(y)][B(y)][A(y)] 


B-1 




(B-1) 


2, and 3, 


(B-2) 


(B-3) 


(B-4) 


(B-5) 



